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1 rROSS-REFERF^rF TO RELAT T?r> APPTJCATIONS 
10 Priority is claimed under 35 U.S.C. 119(e) to U.S. provisional application Serial Nos. 

60/411,718; 60/431,882; and 60/479,970 filed on September 18, 2002; December 9, 2002; 
and June 18, 2003, respectively. The contents of these priority applications are hereby 
incorporated by reference in their entireties. 

15 2 TTTCT ,D OF TH E INVENTION 

The invention relates to methods and compositions for the diagnosis and treatment of 
neuropsychiatric disorders, including major depressive disorder (MDD), bipolar affective 
disorder (BAD), schizophrenia and psychotic depression. More specifically, the mvention 
relates to gene signatures, that are associated with electroshock therapy and are usefial for 

20 diagnosing and treating such neuropsychiatric disorders. 

3 RArKGROUND THE INVENTION 

Current treatments for the psychiatric depression, such as that associated with major 
depressive disorder (MDD), bipolar affective disorder (BAD), and psychotic depression, 
25 comprise several classes of antidepressant drugs. Unfortunately, many patients with these 
forms of depression fail to respond to such drug therapies. In such cases, electroconvulsive 
therapy (ECT) may be the only avaiable treatment, hideed, ECT does produce a high level of 

,1 . .1 . trMtmpnt rppimens have failed. 

success in patients wnere oiner piiiuuiawuiug^v-. 
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Psychiatric depression associated with conditions such as major depressive disorder 
(MDD) and psychotic depression are treated with several class of antidepressant drugs, all of 
which require several weeks to be effective. Unfortunately, many people do not respond to 
these drug therapies. In such cases, electroconvulsive therapy (ECT) remains the treatment. 
Although severe, ECT is recognized as a highly effective and rapid treatment, and it produces 
a high degree of success in patients where other therapeutic regimens have failed. 

Receptor interactions and other pharmacological actions of many chemical 
antidepressants have been well studied and are characterized to at least some degree. By 
contrast, the mechanism or mechanisms of ECT remain poorly understood. Molecular 
effects of electroconvulsive seizure (ECS) in the brains of animals have been studied, as a 
model for ECT (for a review, see Fochtmann, Psychopharmacology Bulletin 1994, 30:321- 
444). The effects of ECS are divers, and include increases in levels of neurotransmitters, 
neuropeptides and synaptic remodeling, including sprouting. Id. ECS also afects several 
brain regions, including the hippocampus, frontal cortex, neostriatum, entorhinal cortex, 
temporal-parietal cortex, and monoaminergic nuclei that project to these areas. A summary 
of the effects reported in these regions is provided in a Table, attached as Appendix to this 
specification (see, Section 8.1. infra)- These changes include, inter alia: increases in 
tyrosine kydroxylase in the locus coeruleus and diverse monoamine nerve terminal regions 
(Masserano et al. Science 1981, 214:662-665); decreases in ^-adrenergic receptors, 
typically in the dentate gyrus of the hippocampus and in the frontal cortex (Biegon & IsraeU, 
Eur. J. Pharmacol. 1986, 123:329-334); and increases of neurotrophic factors which can 
positively modulate monoaminergic neurotransmission (Altar et al, J. Neurochem. 1994, 
63:1021-1032; Mzmoxmas et ai, J. Neurosci. 1995,15:7929-7939; Martin-Iverson a/.,/ 
Neurosci. 1994, 14:1262-1270; Siucisk et al., Pharmacol. Biochem. Behav. 1997, 56:131- 
137). Others have reported that the mRNA for brain-derived neurotrophic factor (BDNF) is 
upregulated in seizures (Isackson et al. Neuron 1991, 6:937-948; Nibuya et al., J. Neurosci. 
1995, 15:7539-7547; Rocamora et al, Brain Res. Mol Brain Res. 1992, 13:27-33). 

Yet, the above-mentioned changes as well as those summarized in the Appendix, 
infra, were identified by investigations that examined changes of only very limited numbers 
of mRNA species in any brain region. No large-scale gene analysis following ECS has been 
reported. Hendriksen et al, (Eur. J. Neurosci 2001. 14:1475-1484) have described measured 
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changes in the expression of multiple genes in the hippocampus, but only eight days after 

local seizures induced by electrical stimulation had been terminated by pentobarbital. 

Lukasiuk et al (Eur. J. Neurosci. 2003, 17:271-279) report experiments where hippocampal 

and temporal lobe expression profiling was conducted in animals subject to amygdala 

kindling, but only before they developed spontaneous seizures. 

***** 

The citation and/or discussion of a reference in this section and throughout the 
specification is provided merely to clarify the description of the present invention and is not 
an admission that any such reference is "prior art" to the invention described herein. 

4. STTMMARY THE INVENTION 

The present invention overcomes problems, such as those described supra, by 
providing ECS signature genes - i.e., nucleic acids that are differentially expressed in 
individuals undergoing an ECS treatment or in subjects undergoing electroconvulsive therapy 
(ECT). Preferred ECS signature genes of the invention are set forth in the tables at the 
appendix, infra. See, in particular, in Sections 8.3-8.5 below. 

The invention also provides ECS gene signatures, which represent collections of one 
or more ECS signature genes and, preferably, expression levels of these genes which are 
indicative of ECS or ECT. In additiona, the invention also provides methods and algorithms 
for identifying ECS signature genes and for selecting preferred ECS signature genes, e.g., for 
use in an ECS gene signature. 

ECS genes signatures of the invention, and the ECS signature genes they comprise, 
are useful, e.g., for diagnostic and prognostic methods of identifying people who have or who 
are susceptible to neuropsychiatric disorders, including inter alia major depressive disorder 
(MDD), bipolar affective disorder (BAD) and psychotic depression. Such methods therefore 
are also provided and considered a part of the present invention. The invention additionally 
provides therapeutic methods, which use ECS gene signatures and ECS signature genes to 
treat, ameliorate or prevent such neuropsychiatric disorders. The invention fiirther comprises 
screening methods, including MPHTS® methods, that use ECS gene signatures and/or ECS 
signature genes to identify compounds that are useful in such methods for treating 
neuropsychiatric disorders. 
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The methods of the invention include methods for identifying compounds for treating 
neuropsychiatric disorders, such as schizophrenia, autism, MDD, BAD, schizophrenia and 
psychotic depression. These methods generally comprise steps of contacting a cell or cells 
with a test compound, and determining expression of oneor more signature genes by the cell 
or cells. These signature genes typicaly comprise nucleic acids that hybridize to a nucleic 
acid selected from the group consisting of SEQ ED NOS:l-152 and the complements thereof 
The determined expression of these one or more signature genes is compared to expression of 
those signature genes in a cell or cell that is not contacted with the test compound. Changes 
in expression of the one or more signature genes (compared to expression in cells not 
contacted with the test compound) indicate that the compound is useful for treating the 
neuropsychiatric disorder. 

The invention additionally provides methods for identifying and selecting ECS 
signature genes. In particular, the invention provides methods for selecting one or more 
nucleic acids that are indicative of an effective therapy for treating a neuropsychiatric 
disorder. The method comprises identifying nucleic acids that are differentialy expressed in 
an individual subjected to electroconvulsive seizure (ECS) compared to an individual not 
subjected to ECS. hi various aspects of these methods, the individual may be subjected to 
either actue ECS or chronic ECS, and the nucleic acids identified by this method are 
preferably ones expressed in the brain or in a region of the brain (preferabfy the frontal cortex 
or the hippocampus). In preferred embodiments, a score value is obtained for each of the 
identified nucleic acids. The score value is preferabfy a function of each gene's differential 
expression in mdividuals subjected to ECS and can be determined according to an objective 
method or algorithm, such as the method set forth in Section 8.2 {infra) for selecting ECS 
signature genes that are most informative and therefore most useful for the assays and other 

methods of this invention. 

The invention further provides kits for detecting an ECS gene signature. Kits of the 
invention generally comprise a plurality of oligonucleotides, each of which is capable of 
specifically hybridizing to a different ECS signature gene. For example, in preferred 
embodiments the oUgonucleotides in such kits are capable of specificaUy hybridizng to ECS 
signature genes selected from the genes set forth in the appendix, infra (see, in particular, at 
Sections 8.3-8.5 below), such as SEQ ID NOS:l-152, as well as homologous and 
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complementary sequences thereof. Preferred kits of the invention may comprise, for 
example, oligonucleotide probes that are immobiUzed on a solid surface or support, such as 
in an expression array. In other embodiments, kits of the invention may comprise a plurality 
of oligonucleotide primers, more preferably a plurality of primer pairs, wherein each pair is 
capable of amplifying a particular ECS signature gene (for example in a PCR reaction). Kits 
of the invention may additionally comprise oligonucleotide primers that are capable of 
priming reverse transcription reactions, for generating cDNA from mRNAs of one or more 
ECS signature genes. In still other embodiments, kits of the invention may contain other 
reagents including, for example, nucleotides that are detectably labeled {e.g., for detecting 
the amplification and or hybridization of ECS signature genes), polymerases, and/or buffers. 

Preferred kits of the invention comprise oUgonucleotides that are capable of 
specifically hybridizing to at least five, and more preferably to at least 10, 15, 20. 25, 30, 35, 
40, 45, 50, 60, 70, 80, 90, 100, 1 10, 120, 130, 140, 150 or more signature genes. In other 
embodiments, a kit of the invention may comprise oligonucleotides that specifically 
hybridize to a set of ECS signature genes consisting of 150, 140, 130, 120, 110, 100, 90, 80, 
70, 60, 50, 45, 40, 35, 30, 25, 20, 15, 10, 5 or fewer ECS signature genes of the invention. 



20 



25 



5. RT^TFF DESCRPTTON OF THE DRAWINGS 

Figures lA-lD plot the mean abundances of each of the 8.799 probe sets measured 
with the Affymetrix U34A rodent chip in the frontal cortex or hippocampus of rats exposed 
to either acute or chronic ECS (y-axes) compared to mean abundances of those same probe 
sets in a sham control group (x-axes). Genes denoted by white points on each of the plots are 
ones passed minimum expression restrictions (described in Section 7.1, infra) whose 
expression changed (in ECS rats compared to sham confrols) by either of 1.5 or more or of 
0.67 or less, and having a p values of less than 0.05. Figure 1 A plots gene expression levels 
measured in the frontal cortex of rats exposed to acute ECS compared to sham controls. 
Figure IB plots gene expression levels measured in the frontal cortex of rats exposed to 
chronic ECS compared to sham controls. Figure IC plots gene expression levels measured 
in the hippocampus of rats exposed to acute ECS compared to sham controls, and Figure ID 
30 plots gene expression levels measured in the hippocampus of rats exposed to chronic ECS 
compared to sham controls. 
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Figure 2 shows aplot of a statistical evaluation demonstrating the effect of sample 
size on the number of identified gene changes or "hits" - i.e., measured changes in 
expression with p-values less than 0.01, regardless of the absolute magnitude of change m 
expression level. Groups of 2-10 rats were randomly composed from data of gene expression 
levels in the hippocampus of sham-treated C'sham")and/or acute ECS-treated ("ECS") rats. 
The number of hits was calculated for each microarray analysis of these samples, and the 
process was repeated 100 times for each sample size. Values plotted in Figure 2 are mean 
values ± 1 standard deviation (SD). In the "mixed vs. mixed" condition, samples were 
randomly selected from all acute hippocampus samples, regardless of the freatment. Hie 
"theoretical chance" number of 40 gene changes is based on a p < 0.01 for the 4,000 probes 
detected in the frontal cortex with at least one present call and a mean abundance of > 100. 

Figure 3 shows a Vemi diagram representation of the distribution of changes in 150 
unique probe sets (of which 135 represent unique genes) in the frontal cortex and 
hippocampus after acute and chronic ECS. The total number of genes whose measured 
expression levels changed in each tissue are indicated m parentheses. 

Figures 4A-4D show plots of the measured magnitude in expression levels for genes 
that were differentially expressed in the frontal cortex or hippocampus of rats treated with 
either acute or chronic ECS. Figures 4A and 4B plot changes in expression levels measured 
in the frontal cortex of rats freated with chronic or acute ECS, respectively, compared to 
sham-freated confrols. Figures 4C and 4D plot changes in expression levels measured in the 
hippocampus of rats freated with chronic and acute ECS. respectively, compared to sham- 
treated confrol animals. For each gene whose expression changed (p < 0.05) in freated 
animals compared with sham-freated confrol animals, the ratio of mean expression in the 
ECS group over the control group is plotted in order of the magnitude of changed. Decreases 
in the expression level of genes are plotted below the unity line, whereas increases in 
expression are plotted above the unity line. X-axes in all four plots are aUgned and plotted at 
the same scale so that the width of the graph represents the number of genes affected in each 
direction. 
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6. DETAILED DESCRIPTION 

To date, the identification of therapeutic compounds to treat neuropsychiatric 
disorders has depended almost entirely on serendipity. That is to say, effective drugs and 
5 other therapies for such disorders have traditionally been discovered by chance and not as the 
result of any directed systematic screening method, hideed, the complex polygenetic nature 
of neuropsychiatric disorders, the subtle structural and cellular changes that they entail, and 
the difficulties in diagnosing and monitoring these disorders have made traditional drug 
screening methods extremely difficuU if not impracticable. U.S. patent application 

10 Publication No.2003/0096264 by Altar et al describes screening methods, referred to in that 
appUcation as "Multi-Parameter High Throughput Screening" or "MPHTS," that overcome 
many of these difficulties and are ideally suited for identifying effective and/or promising 
therapeutic compounds to treat neuropsychiatric disorders, including schizophrenia, bipolar 
affective disorder (BAD), autism, major depressive disorder (MDD) and psychotic 

1 5 depression to name a few. 

Briefly, the MPHTS approach pertains to the combination of data generated fi-om 
gene expression profiling coupled with methods for the systematic analysis and/or 
employment of such data. Using the MPHTS methods described herein, large numbers of 
candidate compounds may be screened {e.g., in vitro) to identify ones that are particularly 

20 promising (and, as such, most likely to be suitable) for treating a neuropsychiatric disorder in 
vivo (e.g., in an individual). For descriptive purposes, these assays comprise at least two 
tiers. The first tier involves the determination of genes involved in a particular disorder, 
which is preferably a neuropsychiatric disorder. However, in another embodiment which is 
actually preferred here, the genes identified may be ones that are associated with a particular 

25 model for a neuropsychiatric or other disorder. For instance, the examples infra describe 
experiments that use electroconvulsive seizures (ECS) in rats as a model for 
electroconvulsive therapy. 

The second tier of an MPHTS method involves the implementation of systematic 
methods to screen test compounds. Such screening methods may be either existing assay 

30 platforms that are already known in the art, such as those described in U.S. patent aplication 
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first time. Preferably, however, the screening assays will be automated and/or will be high- 
throughput assays, so that a large number of test compounds can be rapidly screened with a 
minimal amount of labor and effort. 

5 6.1. Definitions 

The terms used in this specification generally have their ordinary meanings in the art, 
withint the context of this invention and in the specific context where each term is used. 
Certain terms are discussed below, or else in the specification, to provide additional guidance 
to the practitioner in describing the compositions and methods of this invention and how they 
10 may be made and used. In accordance with the invention, there may be employed 

conventional techniques in the fields of molecular biology, microbiology, and recombinant 
DNA technology. These techniques are explained more M\y in the literature and are within 
the skill of the art. See, for example, Sambrook, Fitsch & Maniatis, Molecular Cloning: A 
Laboratory Manual, Second Edition (1989) Cold Spring Harbor Laboratory Press, Cold 
15 Spring Harbor, New York (referred to as "Sambrook et al., 1989"); DNA Cloning: A 

Practical Approach, Volumes I and II (D.N. Glover ed. 1985); Oligonucleotide Synthesis 
(M.J. Gait ed. 1984); Nucleic Acid Hybridization (B.D. Hames & S.J. Higgins, eds. 1984); 
Animal Cell Culture (R.I. Freshney, ed. 1986); Immobilized Cells and Enzymes (IRL Press, 
1986); B.E. Perbal, A Practical Guide to Molecular Cloning (1984); F.M. Ausubel et al. 
20 (eds.), Current Protocols in Molecular Biology, John Wiley & Sons, Inc. (1994). 

The term "neuropsychiatric disorder", which may also be referred to as a "major mental 
ilhiess disorder" or "major mental illness", refers to a disorder which may be generally 
characterized by one or more breakdowns in the adaptation process. Such disorders are 
therefore expressed primarily in abnormalities of thought, feeling and/or behavior producing 
25 either distress or impairment of fimction (Le., impairment of mental fimction such as with 
dementia or senility). Currently, individuals may be evaluated for various neuropsychiatric 
disorders using criteria set forth in the most recent version of the American Psychiatric 
Association's Diagnostic and Statistical Manual of Mental Health (DSM-IV). 
Exemplary neuropsychiatric disorders include, but are not limited to, schizophrenia, attention 
30 deficit disorder (ADD), schizoaffective disorder, bipolar affective disorder, unipolar affective 
disorder, and adolescent conduct disorder, 

-8- 

{M:\3235\100m087us2\00058277.DOC MBiiiffliimillimilll } 



f 



As used herein, the term "isolated" means that the referenced material is removed 
from the environment in which it is normally found. Thus, an isolated biological material 
can be free of cellular components; i.e., components of the cells in which the material is 
found or produced. In the case of nucleic acid molecules, an isolated nucleic acid includes a 
5 PGR product, an isolated mRNA, a cDNA, or a restriction fragment. In another embodiment, 
an isolated nucleic acid is preferably excised from the chromosome in which it may be found, 
and more preferably is no longer joined to non-regulatory, non-coding regions, or to other 
genes, located upsfream or downstream of the gene contained by the isolated nucleic acid 
molecule when found in the chromosome, hi yet another embodiment, the isolated nucleic 
10 acid lacks one or more infrons. Isolated nucleic acid molecules include sequences inserted 
into plasmids, cosmids, artificial chromosomes, and the like. Thus, in a specific 
embodiment, a recombinant nucleic acid is an isolated nucleic acid. An isolated protein may 
be associated with other proteins or nucleic acids, or both, with which it associates in the cell, 
or with cellular membranes if it is a membrane-associated protein. An isolated organelle, cell, 
1 5 or tissue is removed from the anatomical site in which it is found in an organism. An 
isolated material may be, but need not be, purified. 

The term "purified" as used herein refers to material that has been isolated under 
conditions that reduce or eliminate the presence of unrelated materials, i.e., contaminants, 
including native materials from which the material is obtained. For example, a purified 
20 protein is preferably substantially free of other proteins or nucleic acids with which it is 
associated in a cell; a purified nucleic acid molecule is preferably substantially free of 
proteins or other unrelated nucleic acid molecules with which it can be found within a cell. 
As used herein, the term "substantially free" is used operationally, in the context of analytical 
testing of the material. Preferably, purified material substantially free of contaminants is at 
25 least 50% pure; more preferably, at least 90% pure, and more preferably still at least 99% 
pure. Purity can be evaluated by chromatography, gel elecfrophoresis, immunoassay, 
composition analysis, biological assay, and other methods known in the art. 

Methods for purification are well-known in the art. For example, nucleic acids can be 
purified by precipitation, chromatography (including preparative solid phase 
30 chromatography, oligonucleotide hybridization, and triple helix chromatography), 

ultracentrifiigation, and other means. Polypeptides and proteins can be purified by various 
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methods including, without limitation, preparative disc-gel electrophoresis, isoelectric 
focusing, HPLC, reversed-phase HPLC, gel filtration, ion exchange and partition 
chromatography, precipitation and salting-out chromatography, extraction, and 
countercurrent distribution. For some purposes, it is preferable to produce the polypeptide in 
a recombinant system in which the protein contains an additional sequence tag that facilitates 
purification, such as, but not limited to, a polyhistidine sequence, or a sequence that 
specifically binds to an antibody, such as FLAG and GST. The polypeptide can then be 
purified from a crude lysate of the host cell by chromatography on an appropriate solid-phase 
matrix. Alternatively, antibodies produced against the protein or against peptides derived 
therefrom can be used as purification reagents. Cells can be purified by various techniques, 
including centrifiigation, matrix separation (eg., nylon wool separation), panning and other 
immunoselection techniques, depletion (e.g., complement depletion of contaminating cells), 
and cell sorting (e.g., fluorescence activated cell sorting or "FACS"). Other purification 
methods are possible. A purified material may contain less than about 50%, preferably less 
than about 75%, and most preferably less than about 90%, of the cellular components with 
which it was originally associated. The "substantially pure" indicates the highest degree of 
purity which can be achieved using conventional purification techniques known in the art. 

A "sample" as used herein refers to a biological material which can be tested, e.g., for 
the presence of one or more polypeptide or nucleic acids. For example, in one embodiment, 
a sample is a sample of nucleic acids from a cell (e.g., mRNA, or nucleic acids derived 
therefrom) and is tested or analyzed for the presence or absence of certain particular nucleic 
acid sequences, corresponding to certain genes that may be expressed by the cell. Such 
samples can be obtained from any source, including tissue, blood and blood cells, including 
circulating hematopoietic stem cells (for possible detection of protein or nucleic acids), plural 
effusions, cerebrospinal fluid (CSF), ascites fluid, and cell cultiire. 

Non-human animals include, without limitation, laboratory animals such as mice, 
rats, rabbits, hamsters, guinea pigs, etc.; domestic animals such as dogs and cats; and, farm 
animals such as sheep, goats, pigs, horses, and cows. A non-human animal of the present 
invention may be a mammaUan or non-mammalian animal; a vertebrate or an invertebrate. 

In preferred embodiments, tiie terms "about" and "approximately" shall generally 
mean an acceptable degree of error for the quantity measured given the nahure or precision of 
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the measurements. Typical, exemplary degrees of error are within 20 percent (%), preferably 
within 10%, and more preferably within 5% of a given value or range of values. 
Alternatively, and particularly in biological systems, the terms "about" and "approximately" 
may mean values that are within an order of magnitude, preferably within 5-fold and more 

5 preferably within 2-fold of a given value. Numerical quantities given herein are approximate 
unless stated otherwise, meaning that the term "about" or "approximately" can be inferred 
when not expressly stated. 

The term "molecule" means any distinct or distinguishable structural unit of matter 
comprising one or more atoms, and includes, for example, polypeptides and polynucleotides. 

10 The term "aberrant" or "abnormal", as appUed herein refers to an activity or feature 

which differs from a normal or activity or feature, or to an activity or feature which is within 
normal variations of a standard value. 

For example, an abnormal activity of a gene or protein refers to an activity which 
differs from the activity of the wild-type or native gene or protein, or which differs from the 

1 5 activity of the gene or protein in a healthy subject. An activity of a gene includes, for 
instance, the transcriptional activity of the gene which may result from, e.g., an aberrant 
promoter activity. Such an abnormal transcriptional activity can result, e.g., from one or 
more mutations in a promoter region, such as in a regulatory element thereof An abnormal 
transcriptional activity can also result from a mutation in a transcription factor involved in 

20 the control of gene expression. 

An activity of a protein can be aberrant because it is stronger than the activity of its 
native coimterpart. Altematively, an activity can be aberrant because it is weaker or absent 
related to the activity of its native coimterpart. An aberrant activity can also be a change in 
an activity. For example an aberrant protein can interact with a different protein relative to 

25 its native counterpart. A cell can have an aberrant activity due to overexpression or 

underexpression of a gene or protein. An aberrant activity can result, e.g., from a mutation 
in the gene, which results, e.g., in lower or higher binding affinity of a ligand or substrate to 
the protein encoded by the mutated gene. 

The term "therapeutically effective dose" refers to that amount of a compound or 

30 compositions that is sufficient to result in a desired activity. 
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The phrase "pharmaceutically acceptable" refers to molecular entities and 
compositions that are physiologically tolerable and do not typically produce an allergic or 
similar untoward reaction (for example, gastric upset, dizziness and the like) when 
administered to an individual. Preferably, and particularly where a pharmaceutical 
composition is used in humans, the term "pharmaceutically acceptable" may mean approved 
by a regulatory agency (for example, the U.S. Food and Drug Agency) or listed in a generally 
recognized pharmacopeia for use in animals (for example, the U.S. Pharmacopeia). 

The term "carrier" refers to a diluent, adjuvant, excipient, or vehicle with which a 
compound is administered. Sterile water or aqueous saline solutions and aqueous dextrose 
and glycerol solutions are preferably employed as carriers, particularly for injectable 
solutions. Exemplary suitable pharmaceutical carriers are described in "Reminington's 
Pharmaceutical Sciences" by E.W. Martin. 

The term "polymer" means any substance or compound that is composed of two or more 
building blocks ('mers') that are repetitively linked together. For example, a "dimer" is a 
compound in which two building blocks have been joined together; a "trimer" is a 
compound in which three building blocks have been joined together; etc. 

The term "polynucleotide" or "nucleic acid molecule" as used herein refers to a 
polymeric molecule having a backbone that supports bases capable of hydrogen bonding to 
typical polynucleotides, wherein the polymer backbone presents the bases in a manner to 
permit such hydrogen bonding in a specific fashion between the polymeric molecule and a 
typical polynucleotide (e.g., single-stranded DNA). Such bases are typically inosine, 
adenosine, guanosine, cytosine, uracil and thymidine. Polymeric molecules include "double 
stranded" and "single stranded" DNA and RNA, as well as backbone modifications thereof 
(for example, methylphosphonate linkages). 

Thus, a "polynucleotide" or "nucleic acid" sequence is a series of nucleotide bases 
(also called "nucleotides"), generally in DNA and RNA, and means any chain of two or more 
nucleotides. A nucleotide sequence frequently carries genetic information, including the 
infomation used by cellular machinery to make proteins and enzymes. The terms include 
genomic DNA, cDNA, RNA, any synthetic and genetically manipulated polynucleotide, and 
both sense and antisense polynucleotides. This includes single- and double-stranded 
molecules; i.e., DNA-DNA, DNA-RNA, and RNA-RNA hybrids as well as "protein nucleic 
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acids" (PNA) formed by conjugating bases to an amino acid backbone. This also includes 
nucleic acids containing modified bases, for example, thio-uracil, thio-guanine and fluoro- 
uracil. 

The polynucleotides herein may be flanked by natural regulatory sequences, or may 
5 be associated with heterologous sequences, including promoters, enhancers, response 

elements, signal sequences, polyadenylation sequences, introns, 5'- and 3'-non-coding regions 
and the like. The nucleic acids may also be modified by many means known in the art. Non- 
limiting examples of such modifications include methylation, "caps", substitution of one or 
more of the naturally occurring nucleotides wifli an analog, and intemucleotide modifications 
10 such as, for example, those with uncharged linkages (eg., methyl phosphonates, 

phosphotriesters, phosphoroamidates, carbamates, etc.) and with charged linkages (e.g., 
phosphorothioates, phosphorodithioates, etc.). Polynucleotides may contain one or more 
additional covalently Unked moieties, such as proteins {e.g., nucleases, toxins, antibodies, 
signal peptides, poly-L-lysine, etc.), intercalators (eg., acridine, psoralen, etc.), chelators 
15 (eg, metals, radioactive metals, iron, oxidative metals, etc.) and alkylators to name a few. 
The polynucleotides may be derivatized by formation of a methyl or ethyl phosphotriester or 
an alkyl phosphoramidite linkage. Furthermore, the polynucleotides herein may also be 
modified with a label capable of providing a detectable signal, either directly or indirectly. 
Exemplary labels include radioisotopes, fluorescent molecules, biotin and the like. Other 
20 non-limiting examples of modification which may be made are provided, below, in the 
description of the present invention. 

Specific non-limiting examples of synthetic nucleic acids envisioned for this 
mvention include, in addition to the nucleic acid moieties described above, nucleic acids that 
contain phosphorothioates, phosphotriesters, methyl phosphonates, short chain alkyl, or 
25 cycloalkyl intersugar linkages or short chain heteroatomic or heterocyclic intersugar 
linkages. Most preferred are those with CH2-NH-O-CH2, CHrN(CH3)-0-CH2, CH^-O- 
N(CH3)-CH2, CH2-N(CH3)-N(CH3)-CH2 and 0-N(CH3)-CHrCH2 backbones (where 
phosphodiester is O-PO2-O-CH2). US Patent No. 5,677,437 describes heteroaromatic nucleic 
acid linkages. Nitrogen linkers or groups containing nitrogen can also be used to prepare 
30 nucleic acid mimics (U.S. Patents Nos. 5,792,844 and 5,783,682). US Patent No. 5,637,684 
describes phosphoramidate and phosphorothioamidate oligomeric compounds. Also 
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envisioned arc nucleic acids having morphotoo backbone structures (U.S. Pat. No. 
5 034 506) in oa>er embodiments, such as the peptide-nucleic acid (PNA) backbone, the 
phosphodiester backbone of the nucleic acid may be replaced with a polyamide backbone, 
the bases being bound directly or ind.rectly to the aza nitrogen atoms of the polyam.de 
5 backbone (Nielsen e, aU Science 254:1497. 1991). Other synthetic nucleic actds may 
contain substituted sugar moiefies comprising one of the following at the X position: OH. 
SH SCH,. F. OCN, 0(Cft)„NH, or 0(CH,).Cft where n is torn 1 to about 10; C, to C. 
lower alkyl. substituted lower alkyl. alkaryl or aralkyl; CI; Br; CN; CF=; OCF,; 0-; S-. or N- 
alkyl O- S- orN-alkenyl;SOCH,;SO.CH,;ONO.;NO^;N,;NHa;heterocycloalkyl; 
10 heter^ycloalkaryl; aminoalkylamino; polyalkylamino; substituted silyl; a fluorescem 

moiety, an RNA cleaving group; a reporter poup; an intercalator, a group for improvmg the 
pharmacokineHc properties of a nucleic acid; or a group for improvmg the pharmacodynatmc 
properties of an nucleic acid, and other substituents having similar properties. Nucletc actds 
may also have sugar mimetics such as cyclobutyls or other carbocyclics in place of the 
15 pentofuranosyl group. Nucleotide units having nucleosides other than adenosine, cyudtne, 
guanosine. thymidine and uridine, such as inosine. may be used in an oUgonucleottde 
molecule. 

The term "oligonucleotMe" refers to a t^ucleic acid, generally of at least 10. preferably 
at least 15, and more preferably at least 20 nucleotides, preferably no more than 100 

20 nucleotides, that is hybridizable to a genomic DNA molecule, a cDNA molecule, or an 
molecule encoding a gene, mRNA, cDNA, or other nucleic acid of mterest. 
OUgonucleotides can be labeled, e.g. , with -P-nucleotides or nucleotides to which a label, 
such as biotin or a fluorescent dye (for example, Cy3 or Cy5) has been covalently 
conjugated. In one embodiment, a labeled oligonucleotide can be used as a probe to detect 

25 thepresenceofanucleicacid. M another embodiment, oligonucleotides (one or both of 
which may be labeled) can be used as PCR primers, either for cloning full length or a 
fragment of a gene, or to detect the presence of nucleic acids encoding a particular gene 
product (eg., to detect the presence of a particular mRNA). hi a further embodiment, an 
oligonucleotide of the invention can form a triple helix. Generally, oligonucleotides are 

30 prepared synthetically, preferably on a nucleic acid synthesizer. Accordingly, 
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oligonucleotides can be prepared with non-naturally occurring phosphoester analog bonds. 

such as thioester bonds, etc. 

A "polypeptide" is a chain of chemical building blocks called amino acids that are 
linked together by chemical bonds called "peptide bonds". Hie term "protein" refers to 
polypeptides that contain the amino acid residues encoded by a gene or by a nucleic acid 
molecule (e.g., an mRNA or a cDNA) transcribed from that gene either directly or indirectly. 
Optionally, a protein may lack certain amino acid residues that are encoded by a gene or by 
an mRNA. For example, a gene or mRNA molecule may encode a sequence of amino acid 
residues on the N-terminus of a protein (i.e., a signal sequence) that is cleaved from, and 
therefore may not be part of, the final protein. A protein or polypeptide, including an 
enzyme, may be a "native" or "wild-type", meaning that it occurs in nature; or it may be a 
"mutant", "variant" or "modified", meaning that it has been made, altered, derived, or is in 
some way different or changed from a native protein or from another mutant. 

A "Ugand" is, broadly speaking, any molecule that binds to another molecule, hi 
preferred embodiments, the Ugand is either a soluble molecule or the smaller of the two 
molecule or both. The other molecule is referred to as a "receptor". In preferred 
embodiments, both a Ugand and its receptor are molecules (preferably proteins or 
polypeptides) produced by cells. Preferably, a Ugand is a soluble molecule and the receptor 
is an integral membrane protein {i.e., a protein expressed on the surface of a cell). The 
binding of a Ugand to its receptor is frequently a step of signal transduction within a cell. 
Exemplary ligand-receptor interactions include, but are not limited to, binding of a hormone 
to a hormone receptor (for example, the binding of estrogen to the estrogen receptor) and the 
binding of a neurotransmitter to a receptor on the surface of a neuron. 

"AmpUfication" of a polynucleotide, as used herein, denotes the use of polymerase 
chain reaction (PGR) to increase the concentration of a particular DNA sequence within a 
mixture of DNA sequences. For a description of PGR see Saiki et al. Science 1988, 
239:487. 

"Chemical sequencing" of DNA denotes methods such as that of Maxam and Gilbert 
(Maxam-Gilbert sequencing; see Maxam & Gilbert, Proc. Natl. Acad. Sci. U.S.A. 1977, 
74:560), in which DNA is cleaved using individual base-specific reactions. 
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■■Enzymatic sequencing" of DNA denotes methods such as that of Sanger (Sanger « 
al. Proc Na,l.Acad.Sci. USA. 1977, 74:5463) and variation, thereofwell known m the art, 
in "a single-stranded DNA is copied and randomly terminated using DNA polymerase. 

A "gene" is a sequence of nucleotides which code for a functional "gene product". 
5 QeneraUy,age„eproductisaf»nctionalprotein. However, a gene product can also be 
another type of molecule in a cell, such as an RNA (e.g., a tRNA or a rRNA). For the 
purposes of the present invention, a gene product also refers ,„ an mSNA sequence whrch 
may be found in a cell. For example, measuring gene expression levels according to the 
invention may correspond to measuring mRNA levels. A gene may also comprise regulatory 
,0 (,• . . non-coding) sequences as well as coding sequences. Exemplary regulatory sequence 
include promoter sequences, which determine, for example, the conditions under whrch the 
gene is expressed. The transcribed region of the gene may also include untrmslated repons 
including mtrons, a 5'-untranslated region (5'-UTR) and a 3^-un,ranslated region (3'-UTR). 
A "coding sequence" or a sequence "encoding" an expression product, such as a 
15 RNA polypeptide, protein or enzyme, is a nucleotide sequence that, when expressed, results 

in the' production of that RNA. polypeptide, protdn or enzyme; f.e.. toe nucleotide sequence 
"encodes" that RNA or it encodes the amino acid sequence for that polypeptide, protem or 

enzyme. 

A ■■promoter sequence" is a DNA regulatory region capable of binding RNA 
20 polymerase in a ceU and initiating transcription of a downsnream (3- direction) coding 

sequence. For purposes of defining the present invention, the promoter sequence ,s bounded 
a, its 3' terminus by the transcription initiation site and extends upstream (5^ dtrectton) to 
include the minimum number of bases or elements necessary to initiate n^nscription at levels 
detectable above background. Within the promoter sequence will be found a bran^rrptron 
25 initiation site (conveniently found, for example, by mapping with nuclease SI), as well as 
protein binding domains (consensus sequences) responsible for the hindmg of RNA 
polymerase. 

A coding sequence is "under the control of or is "operatively associated with 
transcriptional and translational control sequences in a cell when RNA polymerase 
30 transcribes tj,e coding sequence into RNA, which is then trans-RNA spliced (if it contams 
introns) and, if the sequence encodes a protein, is translated into that protein. 
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The term "express" and "expression" means allowing or causing the information in a 
gene or DNA sequence to become manifest, for example producing RNA (such as rRNA or 
mRNA) or a protein by activating the cellular fiinctions involved in transcription and 
translation of a corresponding gene or DNA sequence. A DNA sequence is expressed by a 
cell to form an "expression product" such as an RNA {e.g., a mRNA or a rRNA) or a protein. 
The expression product itself, e.g., the resulting RNA or protein, may also said to be 

"expressed" by the cell. 

The term "heterologous" refers to a combination of elements not naturally occurring. 
For example, the present invention includes chimeric RNA molecules that comprise an rRNA 
sequence and a heterologous RNA sequence which is not part of the rRNA sequence. In this 
context, the heterologous RNA sequence refers to an RNA sequence that is not naturally 
located within the ribosomal RNA sequence. Alternatively, the heterologous RNA sequence 
may be naturally located within the ribosomal RNA sequence, but is found at a location in 
the rRNA sequence where it does not naturally occur. As another example, heterologous 
DNA refers to DNA that is not naturally located in the cell, or in a chromosomal site of the 
cell. Preferably, heterologous DNA includes a gene foreign to the cell. A heterologous 
expression regulatory element is a regulatory element operatively associated with a different 
gene that the one it is operatively associated with in nature. 

The terms "mutant" and "mutation" mean any detectable change in genetic material, 
20 e.g. , DNA, or any process, mechanism or result of such a change. This includes gene 

mutations, in which the structure {e.g., DNA sequence) of a gene is altered, any gene or DNA 
arising from any mutation process, and any expression product {eg., RNA, protein or 
enzyme) expressed by a modified gene or DNA sequence. Tlie term "variant" may also be 
used to indicate a modified or altered gene, DNA sequence, RNA, enzyme, cell, etc.; i.e., 
any kind of mutant. For example, the present invention relates to altered or "chimeric" RNA 
molecules that comprise an rRNA sequence that is altered by inserting a heterologous RNA 
sequence that is not naturally part of that sequence or is not naturally located at the position 
of that rRNA sequence. Such chimeric RNA sequences, as well as DNA and genes that 
encode them, are also referred to herein as "mutant" sequences. 



15 



25 
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"Sequence-conservative variants" of a polynucleotide sequence are those in which a 
change of one or more nucleotides in a given codon position results in no alteration in the 

amino acid encoded at that position. 

"Function-conservative variants" of a polypeptide or polynucleotide are those in 
which a given amino acid residue in the polypeptide, or the amino acid residue encoded by a 
codon of the polynucleotide, has been changed or altered without altering the overall 
conformation and function of the polypeptide. For example, function-conservative variants 
may include, but are not limited to, replacement of an amino acid with one havmg similar 
properties (for example, polarity, hydrogen bonding potential, acidic, basic, hydrophobic, 
aromatic and the like). Amino acid residues with similar properties are well known m the art. 
For example, the amino acid residues arginine, histidine and lysine are hydrophilic, basic 
amino acid residues and may therefore be interchangeable. Similar, the amino acid residue 
isoleucine. which is a hydrophobic amino acid residue, may be replaced with leucine, 
methionine or valine. Such changes are expected to have little or no effect on the apparent 
molecular weight or isoelectric point of the polypeptide. Amino acid residues other than 
those indicated as conserved may also differ in a protein or enzyme so that the percent 
protein or amino acid sequence similarity (e.g., percent identity or homology) between any 
two proteins of similar function may vary and may be, for example, from 70% to 99% as 
determined according to an alignment scheme such as the Cluster Method, wherein similarity 
is based on the MEGALIGN algorithm. "Function-conservative variants" of a given 
polypeptide also include polypeptides that have at least 60% amino acid sequence identity to 
the given polypeptide as determined, e.g., by the BLAST or FASTA algorithms. Preferably, 
function-conservative variants of a given polypeptide have at least 75%, more preferably at 
least 85% and still more preferably at least 90% amino acid sequence identity to the given 
polypeptide and, preferably, also have the same or substantially similar properties of 
molecular weight and/or isoelectric point) or functions (e.g., biological functions or 
activities) as the native or parent polypeptide to which it is compared. 

The term "homologous", in all its grammatical forms and spelling variations, refers to 
the relationship between two proteins that possess a "common evolutionary origin", including 
30 proteins from superfamilies {e.g., the immunoglobulin superfamily) in the same species of 
organism, as well as homologous proteins from different species of organism (for example, 

-18- 

{M:\3235\100m087us2\00058277.DOC iBIIiiillllllllilMiMnB } 



25 



™y„smUgh.chatapolypep«de,e<c.; see.R=eok .Mi.. 0.1.1987.50:667). such prote,,. 
(and tor encoding nucleic acids) have sequence homology, as reflected by d,eir sequence 
similarity, whether in terms of percent idenmy or by the pre^nce of specific restdues or 

motifs and conserved positions. 
5 The term "sequence similarity", in all its grammatical fomts, refers to the degree of 

identity or correspondence between nucleic acid or amino acid sequences that may or may 
not share a common evolutionary origm (see. Reeck supra). However, in common 
usage and ii> the instant application, the tem. "homologous", when modified with an adverb 
such as "highly", may refer to sequence similarity and may or may not relate to a common 

10 evolutionary origin. 

In specific embodiments, two nucleic acid sequences are "subst^tially homologous 
or "substantially similar" when at least about 80o/„, and more preferably at least about or 
at least about 95% of the nucleotides match over a defined length of the nucleic acid 
sequences, as determined by a sequence comparison algorithm known such as BLAST, 
15 FASTA, DNA Strider, CLUSTAL, etc. An example of such a sequence is an allelic or 
species variant of the specific genes of the present invention. Sequences that are 
substantially homologous may also be identified by hybridization, e.g., in a Southern 
hybridization experiment under, e.g. , stringent conditions as defined for that particular 
system. 

20 Similarly, in particular embodiments of the invention, two amino acid sequences are 

"substantially homologous" or "substantially similar" when greater than 8O0/0 of the amino 
acid residues are identical, or when greater than about 90% of the amino acid residues are 
similar (i.e., are fimctionally identical). Preferably the similar or homologous polypeptide 
sequences are identified by aUgmnent using, for example, the GCG (Genetics Computer 
25 Group, Program Manual for the GCG Package, Version 7, Madison Wisconsin) pUeup 
program, or using any of the programs and algorithms described above {e.g., BLAST, 

FASTA, CLUSTAL, etc.). 

The terms "array" and "microarray" are used interchangeably and refer generally to 
any ordered arrangement (e.g., on a surface or substrate) or different molecules, referred to 
30 herein as "probes". Each different probe of an arrays specifically recognizes and/or binds to 
a particular molecule, which is referred to herein as its "target". Microarrays are therefore 
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useful for simultaneously detecting the presence or absence of a plurality of different target 
molecules, e.g., in a sample. In preferred embodiments, arrays used in the present invention 
are "addressable arrays" where each different probe is associated with a particular "address". 
For example, in preferred embodiments where the probes are immobilized on a surface or a 
5 substrate, each different probe of the addressable array may be immobilized at a particular, 
known location on the surface or substrate. The presence or absence of that probe's target 
molecule in a sample may therefore be readily determined by simply determining whether a 
target has bound to that particular location on the surface or substrate. 

In various embodiments, an array of the invention may comprise a plurality of 
10 different antibodies that each bind to a particular target protein or antigen. More preferably, 
however, the methods of the invention are practiced using nucleic acid arrays (also referred 
to herein as "transcript arrays" or "hybridization arrays") that comprise a plurality of nucleic 
acid probes immobilized on a surface or substrate. The different nucleic acid probes are 
complementary to, and therefore hybridize, to different target nucleic acid molecules, e.g., m 
15 a sample. Thus such probes may be used to simultaneously detect the presence and/or 
abundance of a plurality of different nucleic acid molecules in a sample, including the 
expression of a plurality of different genes; e.g., the presence and/or abundance of different 
mRNA molecules, or of nucleic acid molecules derived therefrom (for example, cDNA or 
cRNA). 

20 A nucleic acid molecule is "hybridizable" to another nucleic acid molecule, such as a 

cDNA, genomic DNA, or RNA, when a single stranded form of the nucleic acid molecule 
can amieal to the other nucleic acid molecule under the appropriate conditions of temperature 
and solution ionic strength {see Sambrook et al, supra). The conditions of temperature and 
ionic strength determine the "stringency" of the hybridization. For preliminary screening for 

25 homologous nucleic acids, low stringency hybridization conditions (e.g., 5x SSC, 0.1% SDS, 
and no formamide; or 30o/o formamide, 5x SSC, O.S^/o SDS) may be used. Alternatively, 
hybridizations may also be performed under conditions that are relatively more stringent, 
such as moderately stringent hybridization conditions {e.g., 40% formamide, with 5x or 6x 
SCC) or high stringency hybridization conditions (e.g., 50% formamide, 5x or 6x SCC). 

30 SCC is a buffer solution commonly used for nucleic acid hybridizations and comprises 
0.15 M NaCl, 0.015 M Na-citrate. 
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Hybridization requires tl«t the two nucleic acids contain complementary sequences, 
although dependingon the stringency ofthehybridization,mismatcl,eabetwe«,bases are 

possible. The appropriate stringency for hybridizing nucleic »:ids depends on the len^ot 
Ure nucleic ac.ds and the degree of complementation, variables well known m the art. Tire 
5 greater the degree of similarity or homology between two nucleotide sequences, the ^eater 
the value ofT^forhybridsofnucleic acids havingthose sequences. TherelativestabtUty 

(corresponding to higher of nucleic acid hybridizations decreases in the following on^er 
RNA-RNA, DNA:RNA. DNA:DNA. For hybrids of greater than 100 nucleotides m length, 
equarions for calculatingT^havebeen derived (...Sambrook.ra(.,.»pra,9.50-9.51> For 

,0 hybridizationwithshorternucleicacids.i.e.,oligonucle„tides.,hepositionofmis.^tches 

becomes more important, and the length of the oli^nucleotide determines its spectoty 
Sambrook .,a;..»,ra, 11.7-11.8). Aminimum length forahybridizablenucleicacrdts a. 

leas. aboutlOn„cleo.ides;preferab.yat least abou.l5nucleorides;andmorepreferab.y.he 

length is at least about 20 nucleotides. 

15 Suitable hybridization conditions for oligonucleotides (e.g.. for oligonucleonde 

probes or primers) are typically somewhat differem than for Ml-lengU. nucleic acds (e.g 
full-length cDNA). because of the oligonucleotides' lower meMng temperature. Because the 
melting temperature of oligonucleorides wiU depend on the length of the oligonucleotrde 
sequences involved, suitable hybridization temperatures will vary depending upon the 

20 oligon«cleotidemoleculesused.Exemplaryteraperaturesmaybe37X(forl4-base 

oUgonucleotides), 48 «C (for 17-base oligonucleotides), 55 "C (for 20-base ohgonucleottd^) 
and 60 X (for 23-base oligonucleotides). Exemplary suitable hybridization condihons for 
oligonucleotides include washing in 6x SSOO.05% sodium pyrophosphate, or other 
conditions that afford equivalent levels of hybridization. 
25 Preferably, nucleic acid molecules m the present invention are detected by 

hybridization to probesofamicn^array. Hybridization and washconditions are the^fore 

preferably chosen so that the probe "specifically binds" or "specifically hybridizes" to a 
specific t^get nucleic acid, h, other words, the nucleic acid probe preferably hybridtzes, 

duplexes or binds to a target nucleic acid molecules having a complementary nucleotide 
30 sequence, but does not hybridize to a nucleic acid molecules having a non-complementary 
sequence. As used herein, one polynucleotide sequ^tce is considered complcucntary to 
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anott-er when, if toe shorter otfte polynucleotides is less to or equal to abou, 25 bases, 
there are no mismatches nsingsundara base-pairing rules. Itthe shorter of are two 

polynucleotides is longer than about 25 bases, there is preferably no more than a 5% 
mismatch. Preferably, the two polynucleotides are perfectiy complementary (i.e., no 
5 mismatches). In can be easily demonstrated that particular hybridization condtUons are 
suitable for specific hybridization by carrying out the assay using negative controls. See, for 
example. Shalon et a/., Geno,^ ResearcH 1996. 639-645; and Chee e, ai. Science 1996. 

274:610-614. , ^ , ^. 

Optimal hybridization conditions for use with mioroarrays wiU depend on the length 
,0 (eg. „ligonucleotidever^polynucleotidegreaterthanabo„t200bases)andtype(..«., 
RNA DNA,PNA,e<c.)ofprobeandtarge,nucleicacid. General parameters for spectfic 
<ie Stringent) hybridization conditions are de^bed above. For cDNA mioroarrays, such 
as ^ose described by Schena.r<,/. iFroc. Na,l. Aca,. Sc. USA 1996. 93:.<XS.4) typica. 
hybridization conditions comprise hybridizing in 5x SSC and 0.2% SDS at 65 X for about 
,5 four hours, followed by washes at 25 «C in a low sttingency wash buffer (for example, Ix 
SSC and 0 2% SDS), and about 10 minutes washing at 25 "C in a high sttingency wash 
buffer (for example, 0. Ix SSC and 0.2% SDS). UseftU hybridization conditions are also 
provided, e.g.. in Tijessen, Hyl,riclm,ion .UH Nucleic Acid Probes. Elsevier Scences 
Pubhshers (1 996), and Kricka, Nonisoiopic DNA Probe Techniques. Academic Press. San 

20 Diego OA (1992). . 

The term "expression profile" or "gene signature" refer, generally, to any descnption 
or measurement of the genes and/or nucleic acids that are expressed by a cell or organrsm 

nnder particular conditions. For example, an expression profile may be measured under 
particular conditions of ^wtt.. for example ataparticulartemperanne, in tire presence or 

25 absence of particular growti, media, and/or in the pres^ce or absence of particular numents. 
h, preferred embodiments, gene signatures may be obtained, e.g.. for cells or tissues that are 
derived from an individual or individuals having a neuropsychiatric disorder. Gene 
signati^es may also be obtained for a cell or organism exposed to one or mote particular 
drugs or other compounds, such as for acell or organism exposed to a known therapeutic 

30 compound(eg..witt.alcnownusefortteatiuganeu,»psychia.ricdisorder)orforacellor 

organism exposed to a "test" or "candidate" compound ^e.g., as part ox . .. .... u...... ^ 
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expression profile or gene signature may comprise a description of particular genes that are 
expressed by a cell or organism, a description of the level or abundance with which genes are 
expressed in a cell or organism, or both. Accordmgly, the term "signature.gene" is used 
herein to refer to a gene that may be used, either alone or with other genes (e.g., as part of a 
gene signature) to characterize a particular condition such as the presence or absence of a 

neuropsychiatric disorder. 

Preferably, an expression profile will comprise a list of different mRNA species that 
are expressed by a cell and their relative abundances. For example, mRNA abundances can 
be measured using a microarray, as described below. In more preferable embodiments, 
nucleic acids (e.g., mRNA) expressed by a cell are reversed transcribed into either cDNA or 
cRNA, and the abundances of the cDNA and/or cRNA molecules are measured. 

6.2. Multi-Parameter High Thron ^^bput Screening (MPHTS) 
MPHTS methods are described in U.S. patent application Publication No. 
2003/0096264 by AUar et al. Such methods generally comprise the following five elements. 
The skilled artisan will appreciate, however, that the invention may be practiced omitting one 
or more of these elements and without executing the recited elements in any particular order. 
For example, in certain embodiments, some of the below-described elements may be 
obtained from another source, such as from an online database. The invention may therefore 
be practiced without necessarily performing each of these elements, e.g. , as a separate step in 
a screening method. 

First, gene-signatures are obtained or provided by measuring expression levels for a 
plurality of genes in cells or tissues derived from an individual having a neuropsychiatric 
disorder or from a suitable model. For instance, in the Examples, infra, electroconvulsive 
seizures (ECS) induced in rats are used as models for electroconvulsive therapy (ECT) in 
human patients. Gene signatures are therefore obtained by measuring expression levels for a 
plurality of genes in the cells or tissues (more specifically in frontal cortex and hippocampus 
tissues) of animals subjected to ECS. hi alternative embodiment, the cells and/or tissues are 
brain cells or tissues can be derived from human psychiatiic patients (for example, in post 
mortem tissue samples). However, brain and other neuronal cells or tissues from other 
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species of organisms may also be used, such as from a mouse, a rat, a primate or another 
speciesofmammal.Preferably,the organism from which the brain cells or tissue are denved 

represents an acceptable animal model for a neuropsychiatric disorder. Preferably, the 
expression levels measured in the cells or tissues are compared to expression levels from 
nomial cells or tissues e. , brain cells or tissues from healthy individuals, not affected by a 
neuropsychiafric disorder) to identify particular genes that are differentially expressed m 
ceUs from an individual having a neuropsychiatric disorder compared to one who does not 
have a neuropsychiatric disorder. 

Second, gene-signatures may also be obtained or provided by measuring expression 
levels for a plurality of genes in cultured neuronal cells or tissues {e.g. , in cultured neurons 
that are derived from neural stem cells or from other neuronal cell lines). Human neurons 
and/or neuronal cell lines are particularly preferred. However, the cells may be obtained or 
derived from any species of organism, particularly a mammalian species such as a mouse, a 
rat or a primate. Similarly, the cultured neuronal tissues may also be obtained from any 
species of mammal, such as from a rat, a mouse, a primate or a human. 

For example, and not by way of limitation, a mouse neuroblastoma cell line may be 
used in such methods. Such cells are readily available, ..g., from the American Type Culture 
Collection ("ATCC", Manasas Virginia). See, for example, ATCC Accession No. CRL- 
2263 As another non-limiting example, U.S. provisional patent application serial no. 
60/299 066 filed on June 18, 2001 describes the use of rat neuronal cell cultures to evaluate 
neuropsychiatric drugs. Such cells may also be used in the MPHTS methods of this 
invention. 

Third, drug signatures may also be obtained or provided by measuring expression 
levels for a plurality of genes in cultured neuronal cells or tissues that are freated with a 
therapeutic compound. The cultured cells may be any type of neuronal cell or cell hues 
described supra for obtaining gene-signatures from a cell line. Similarly, any of the types of 
tissue cultures described, supra, may also be used to obtain drug signatures. Preferably, the 
) drug signatures are signatures for compounds that are known to be effective for treatmg a 
neuropsychiafric disorder. Exemplary compounds may mciude valproate, buspironc, lithrom, 
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carbamazepine, clozapine, olanzapine, haloperidol, secretin and vasoactive intestinal 
polypeptide (VIP), to name a few. Exemplary drug signatures, which were obtained from 
broth rat and human neuronal cells treated with therapeutic compounds, are provided in the 
Examples, infra. Other drug signatures may be readily obtained by those skilled in the art. 

Fourth, expression levels for the plurality of genes are obtained or provided in 
neuronal cells that are contacted with a test compound (referred to here as a "drug 
candidate"), and these expression levels may then be compared to expression levels from 
gene signatures obtained for the neuropsychiatric disorder (as described in the first element, 
supra) and/or to drug-signatures obtained the known therapeutic compound (as described in 
the third element, supra). In preferred embodiments, expression levels or "signatures- 
obtained from a test compound are also compared to expression levels when the cell or cell 
line is not contacted with the test compound or any other drug (described in the second 
element, supra). 

Generally speaking, the "signature" or expression levels obtained when the neuronal 
cells are contacted with a test compound are compared to the gene signatures of the cells 
when they are not contacted with any test or therapeutic compound (i.e., the gene signature 
obtained as element two, described supra) to identify changes in the expression level(s) for 
particular genes. Similarly, the drug-signature (obtained as described, supra, for element 
three) is also compared to the neuronal cell lines gene signature, to identify particular genes 
whose expression levels change when the cells are contacted with the therapeutic compound. 
In instances where changes in expression levels when the cells are contacted with the test 
compound are identical (or at least similar) to changes in expression levels when the cell are 
contacted with the known therapeutic compound, then the test compound is identified as a 
candidate compound for treating the neuropsychiatric disorder. Thus, using these screening 
methods a skilled artisan is able to rapidly and inexpensively identify compounds that are 
most promising as novel neuropsychiatric drugs, while eliminating compounds that show 
little promise and/or are unlikely candidates for treating a neuropsychiatric disorder. 

In preferred embodiments of the invention, changes in expression levels when the 
cells are contacted with the test compound may also be compared to gene signatures obtained 
for the particular neuropsychiatric disorder of interest {i.e., to the gene signatures obtained as 
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desoribed,™pra.forteflrs.el»n«,). Preferably, a test compound tot is idenufied as a 
candidate therapeutic compound will alter the expression of "signature gene" in a way that ts 
opposite or contrary to the expression observed in the disorder's gene signature. For 
example, where a particular gene is expressed at abnormally high levels to cells or tissues 

5 ftom individuals affected by the particular neuropsychiatric disorder (compared to expression 
levels m ceUs or tissues from individuals no. affected by the disorder), a candid^e compound 
identified in ftese screentog mefltods will preferably inhibit that gene's expression (i.e. , the 
gene is preferably expressed at lower levels when the cells are contacted wim the test 
compound, compared to its expression when the cell is not contacted with the test 

10 compound). 

Fifth, as an optional element of the invention, drug candidate or candidate compounds 
that are identified as described, supra, may be further optimized, e.g., to account for 
individual genetic variability. 

As indicated above, the MPHTS assays of the invention are useful as an inexpensive 

and rapid initial screening to quickly identify compounds that are most promising as 
neuropsychiatricdrugs,whilequickly eliminating compounds that show littlepromise and/or 

are unlikely candidates for treating a neuropsychiatric disorder. In preferred embodiments, 
20 the MPHTS assays are used to identify candidate compounds for treating bipolar affective 
disorder (BAD), depression, schizophrenia and autism. However, the assays are by no means 
limited to these particular disorders, and may be readily adapted to identify candidate 
compounds for treating any neuropsychiatric disorder. Other exemplary, preferred 
neuropsychiatric disorders for which these assays may be used mclude anxiety disorders, 
25 eating disorders, addictive disorders and Attention Deficit Hyperactivity Disorder (ADHD). 

Classes of compounds that may be identified by such screening assays mclude, but 
are not limited to, small molecules (eg. , organic or inorganic molecules which are less than 
about 2 kd in molecular weight, are more preferably less than about 1 kd in molecular 
weight and/or are able to cross the blood-brain barrier or gain entry into an appropnate cell, 
as well as macromolecules (e.g. , molecules greater than about 2 kd in molecular weight). In 
preferred embodiments, commercially available compound libraries may be purchased and 
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screened in an MPHTS assay of the invention. Examples of preferred libraries include 
TOCRIS (Tocris Cookson. Ltd. Avonmouth Bristol. United Kingdom), SIGMA RBI (Sigma 
AUdrich Inc., St. Louis MO), ChemBridge (ChemBridge Corp.. San Diego CA), Chemdiv 
(ChemDiv Inc.. San Diego CA) and Prestwick (Prestwick Chemical. Inc.. Washington DC), 

5 to name a few. 

The selection of appropriate small molecule compound concentrations for the 
treatment of cells in vitro or for dosing of animals in vivo is preferred to discriminate 
between physiological and toxicological effects of a given compomid. As an initial means 
for detemiining the deleterious effects of a compound or set of compounds, cells may be 
10 seeded {e.g., in multiple-well plates) and treated with a range of compound concentrations. 
The compounds' effect (e.g., its cytotoxic or apoptotoic effect) may then be gauged, e.g., 
using commercially available kits and routine methods well known in the art. 

Compounds identified by these screening assays may also include peptides and 
polypeptides. For example, soluble peptides, fusion peptides members of combinatorial 
15 libraries (such as ones described by Lam et al. Nature 1991, 354:82-84; and by Houghten et 
al. Nature 1991, 354:84-86); members of libraries derived by combinatorial chemistry, such 
as molecular libraries of D- and/or L-configuration amino acids; phosphopeptides, such as 
members of random or partially degenerate, directed phosphopeptide libraries (see, e.g., 
Songyang et al.. Cell 1993, 72:767-778); antibodies, including but not Ihnited to polyclonal, 
20 monoclonal, humanized, anti-idiotypic, chimeric, or single chain antibodies; antibody 

fragments, including but not limited to FAb, F(ab')2. FAb expression library fragments and 

epitope-binding fragments thereof. 

The compounds used in such screening assays are also preferably essential pure and 
free of contaminants which may, themselves, alter or influence gene expression. Compound 

25 purity may be assessed by any number of means that are routine in the art, such as LC-MS 
and NMR spectroscopy. Libraries of test compomids are also preferably biased by using 
computational selection methods which are routine in the art. Tools for such computational 
selection, such as Pipeline PilotJ (Scitegic Inc., San Diego, California) are commercially 
available. Tlie compounds may be assessed using rules such as the "Lipinski criteria" (see, 

30 Lipinski et al.,Adv. DrugDeliv. Rev. 2001, 46:3-26) and/or any other criteria or metrics 
commonly used in the art. 
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6.3. p^^r^r.tinn of Neiirmni ""'^ T'i^^"^ Samples 



10 



15 



20 



25 



S..i„to„»a™pte.. in certain linuted embodiments, brain cells and tiss»« for us. 
i„ the MPHTS methods of this invention may be obtained from individuals (e.g., fix.m 
patients) in a biopsy. However, those skilled in the art will recognize that brain surgenes 
permitting a biopsy are relatively rare and primarily involve surgical excisions (e.g., for the 
treatmentof epilepsy) rather thanbrainregions relevant to neuropsychiatric disorder such as 

schizophrenia or bipolar affective di^rder. In certain embodiments, however, useM dtsease 
profiles may be obtained from cultured peripheral nervous system neurons, such as 
rhinoneuroepifltelial cells. Such cells may be readily obtained from a nasal biopsy, and 
disease profiles from such cells may be used to identify changes in gene expression that are 
associated with neuropsychiatric disorders such as schizophrenia. 

to preferred embodiments, brain cells or Ussues used in the methods of this invenuon 
are instead obtained post-mortem, e.g., from cadavers of individuals who had or extabtted 
symptoms of a neuropsychiatric disorder during their Ufetime. 

Those skilled in the art wiU readily appreciate dmt a large number of carefully 
collected brain tissue samples should p^ferably be obtained to assure statistical reliabtUty 
(see, for example, Torrey e, al., SchizopHr Res. 2000, 44:151; Bahn e, al.J. Chen,. 
Weur—y 2001, 22:79-94; and Vawter e, ai. Brain Res. Bull 2001, 55:641-«50). Thrs 
is particularly desirable where there is considerable heterogeneity in patient age to penrnt 
accounting for age-associated variables (for example, progressive brain degeneration, wtach 
may also occur in schizophrenia). However, smaller samples may be used, e.g. , for 
preUminary screenmg assays where statistical reliabiUty may not be as essenttal. It ts also 
preferable that the samples be matched, eg.. acconJing to the patients' age, sex, cause of 
death and post-mortem intaval. Tire brain samples used preferably are not acuir«l from 
eadavers under circumstances that might fltemselves affect the quality ofthe cells or ttssues 

acquired For example, samples obtained following a prolonged moribund state, a coma, 
hpoxia, pyrexia or stroke preferably are not used in MPHTS methods of the invention. A 
, skiUed anisan may readily recognize such compromised, ante mortem states, e.g., from the 
extent of brain acidosis. Generally, measured postmortem u..u. v..... 
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abou,6.4tadioa« to. ftetissuehasbeenrtjected to suchacompromised ante mortem s,a,e 
and should not be u«d. In addition, the po^ortem tissue pH value is also critical to the 
integrity of mRNA obtained ftom the tissue. „ , k 

is understood that a reliable psychiatric diagnosis and cause of dea^ should also be 
obtained or determmed for the individual. It is, moreover, additionally preferably to identtfy 
fectors such as concomitant medical conditions, medications taken during the patrenfs 
lifetime (particularly immediatelyprior to death), surgical treatments (includrng cancer 

treatmentsjand substance abuse for eachpatient. The hemisphere and regionof the bratn 

from which each sample is taken is also preferably noted and recorded. 

Generally, samples that have been subject to such conditions as may affect the 
reliabilityofgenecxpressionmeasurementsshouldnotbeused. However, in many 

situations the skilled ariisan will recop^ize that such fectors may be sufficiently controlled 
forandthesample,therefore,acceptableforuseinMPHTS.Insuchcases,however,,fs 

preferable ^d often essential that the samples be appropriately matched. As an example, and 
, not by way of limitation, it is recognized that smokmg alters the expr^sion of many genes m 
the hippocampus,are»o„ofthe brain that is also associated with schizophrenia (Wang 

Ab. Soc.Neurosci.2OOl.2n However, the overlap between genes whose expresston 
leJels have been reported as altered bythose two condirions is beUeved to be minimaKsee. 

Wang ai, supra). Therefore, it may be possible to practice MPHTS methods of the 
0 invenrion using s^ples from smoking or non-smoking individuals, provided the samples are 

appropriately matched. 

Those skilled in the art wiU also appreciate that the levels and quality of RNA 
extracted ftom post-mortem samples may be influenced by factors such as the post mortem 
interval (U., the time interval between death and RNA extraction), the refrigeratron t,me 
,5 (,-. therimeintervalftomdeathtopatientstorageinacoldenvironmentXthestoragetrme 

(i e ' the duration of time during which the cadaver is refrigerated). Accordingly, it is 
preferably that such factors be appropriately controlled and that the steps of RNA extracuon 
ftom these tissue samples be as efHcient as possible. In particularly pr.fe.ed embodiments, 
are brain or tissue samples are unfixed (U. , are not treated with protein cross-linker, such as 
30 formalin)andhavenotbeenthawedmorethanonce. 
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In a preferred embodiment, samples of brain tissue may be obtained, e.g., post- 
mortem from cadavers of individuals who (during their lifetime) suffered from or exhibited 
symptoms of a neuropsychiatric disorder. However, single neurons or groups of 
homogeneous neurons may also be extracted from such cadavers, e.g., by laser capture 
5 microdissection (LCM). Using RNA amplification, gene expression profiles may be 

measured for these single cells as well (see, e.g., Eberwine et al, Proc. Natl. Acad. Sci. 1992, 
89:30130-30134; mdLuo et al., Nature Med. 1999,5:117-119). Expression profiles 
obtained from these cells will therefore be particular for the particular cell types extracted, 
and may ultimately provide gene expression profiles that are more clearly ascribed to the 
10 particular cell population. Such gene profiles will typically be more robust, and therefore 
preferable, for evaluating a drug response. 

Brain cells or tissues obtained from animals may also be used. For example, tissue or 
samples from animal models for a neuropsychiatric disorder may be used to model disease 
profiles for that disorder. Alternatively, expression profiles may be obtained from brain cells 
15 or tissues obtained from animals treated with a known anti-psychotic drug or with a test 
compound. In addition, cells from a fransgenic animals may be employed, in which one or 
more genes relevant to a neuropsychiafric disorder have been altered, over-expressed or 
"knocked-out". High throughput in vitro screening of candidate compounds may then be 
carried out using neuronal cells obtained or derived from such a transgenic animal. 



20 



Neuronal cells. In preferred embodiments, the MPHTS methods of the invention also 
used cultured cells or cell lines to screen for candidate therapeutic compounds. Preferably, 
the cells are ones having an expression profile that is typical of neuronal cells or, 
alternatively, they may be cells which can be manipulated to produce an expression profile 
25 typical of neuronal cells. The cells or cell lines used will also, preferably, give rise to 
reproducible changes in their gene expression profiles when contacted with known 
antipsychiatric drugs (for example, valproate). In a particularly preferred embodiment, these 
changes will be opposite changes that are observed in the disease signature. That is to say, in 
such embodiments, genes (or their homologs) normally expressed at higher levels in the 
30 disease signature are preferably expressed at lower levels in cells or cell lines contacted with 
the known antipsychiatric drug, and vice-versa. 
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to a prefemd embodiment, pluripotent neuronal stem cell lines are used in these 
asp«=ts of the invention. Such cell lines are well known in the art, and methods to induce or 
enhance the differentiation of such stem cell lines have been described. For example, U.S. 
Provtsional Patent Application Publication Nos. 2003/0082802 Al and 2003/00.3192 Al 
5 both describe methods for inducing differentiation in nemonal stem ceUs by exposure to 
chemicals (for example, valproate and buspirone). to other embodiments, such cells may be 
differentiated, e.g., using antisense strategies and/or routine techniques of molecular btology 
to develop stable, transfected cell Unes. Alternatively, however, cells or cell lines may also 
be obtained from patients having a neuropsychiatric disorder of interest. 
10 A skilled artisan will readily appreciate that cells or cell culmres used in the methods 

of this invention should be carefully cootix>Ued for parameters such as the cell passage 
number, cell density (e.g., in microplate wells), the method(s) by which cells are dispensed, 
and growth time after dispensing. It is also preferable to repeat mRNA and/or protem 
expression levels measured for a cell or cell line under particular conditions, to confirm that 
1 5 the measured levels are reproducible. 

64 ^^pai jiiring Gq ie P.xnression 

MPHTS methods generally, including the particular methods of this invention, may 
be implemented using any mefliod suitable for measuring changes in the gene expression of a 

20 cell or cells. Such methods are well known and routinely used in the art. In preferred 

embodiments, methods are used that permit tite simultaneous measurement of expression tor 
a plurality of genes (eg. , for at least 10, more preferably for at least 100, still more preferably 
for at least 150). For example, in particularly prefenred embodiments expression profiles are 
measnrol using "transcript arrays" or -microarrays," as described below. However, any 

25 technique ti>at is capable of measuring gene expr^ion may be used and tire methods of this 
invention are not limited to the use of nucleic acid microarrays. For instance, gene 
expression can also be measured in alternative embodiments by using a reverse transcription 

polymerase chain reaction (RT-PCR). 

Systems and kits for implementing such assays are commercially available from a 
30 number of suppliers, including Affymetrix (Santa Clara, Calfomia), Agilent (Palo AUo. 
California), Promega (Madison, Wisconsin), Xanthon (Research Triangle Park, North 
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Carolina), Illumina (San Diego, California), Chromagen (San Diego, California), Third Wave 
Technologies (Madison, Wisconsin), Aclara (Mountain View, California), Beckton 
Dickinson & Co. (Franklin Lakes, New Jersey) and Luminex (Austin, Texas) to name a few. 

5 Transcript arrays. In a preferred embodiment the present invention makes use of 

"transcript arrays" (also called herein "microarrays"). Transcript arrays can be employed for 
analyzing the steady state level of mRNAs in a cell, and especially for comparing the steady 
state levels between two cells, such as a first cell that has been exposed to a drug, drug 
candidate or other compound, and a second cell that has not been treated. 
10 m one embodiment, transcript arrays are produced by hybridizing detectably labeled 

polynucleotides representing the mRNA transcripts present in a cell (e.g., fluorescently 
labeled cDNA synthesized from total cell mRNA) to a microarray. As explained in the 
definitions, supra, microarray is a surface with an ordered array of binding (e.g., 
hybridization) sites for products of many of the genes in the genome of a cell or organism, 
15 preferably most or ahnost all of the genes. Microarrays can be made in a number of ways, of 
which several are described below. However produced, microarrays share certain 
characteristics. The arrays are preferably reproducible, allowing multiple copies of a given 
array to be produced and easily compared with each other. Preferably the microarrays are 
small, usually smaller than 5 cm^ and they are made fi-om materials that are stable under 
20 binding {e.g., nucleic acid hybridization) conditions. A given binding site or unique set of 
bmding sites in the microarray will specifically bind the product of a single gene in the cell. 
Although there may be more than one physical binding site (hereinafter "site") per specific 
mRNA, for the sake of clarity the discussion below will assume that there is a single site. It 
will be'appreciated that when cDNA complementary to the RNA of a cell is made and 
25 hybridized to a microarray under suitable hybridization conditions, the level of hybridization 
to the site in the array corresponding to any particular gene will refiect the prevalence in the 
cell of mRNA transcribed from that gene. For example, when detectably labeled (e.g., with a 
fluorophore) cDNA complementary to the total cellular mRNA is hybridized to a microarray, 
the site on the array corresponding to a gene {i.e., capable of specifically binding a nucleic 
30 acid product of the gene) that is not transcribed in the cell will have little or no signal, and a 
gene for which the encoded mRNA is prevalent will have a relatively strong signal. 
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In preferred embodiments, cDNAs from two different cells, e.g., a cell exposed to a 
test compound and a cell of the same type not exposed to the compound, are hybridized to 
the binding sites of the microarray. The cDNA derived from each of the two cell types are 
differently labeled so that they can be distinguished. In one embodiment, for example, cDNA 
from a cell treated with a drug is synthesized using a fluorescein-labeled dNTP, and cDNA 
from a second cell, not drug-exposed, is synthesized using a rhodamine-labeled dNTP. When 
the two cDNAs are mixed and hybridized to the microarray, the relative intensity of signal 
from each cDNA set is determined for each site on the array, and any relative difference in 
abundance of a particular mRNA detected. 

In the example described above, the cDNA from the freated cell will fluoresce green 
when the fluorophore is stimulated and the cDNA from the untreated cell will fluoresce red. 
As a result, when the compound has no effect, either directly or indirectly, on the relative 
abundance of a particular mRNA in a cell, the mRNA will be equally prevalent in both cells 
and, upon reverse transcription, red-labeled and green-labeled cDNA will be equally 
prevalent. When hybridized to the microarray, the binding site(s) for that species of RNA 
will emit wavelengths characteristic of both fluorophores. hi contrast, when the cell is 
exposed to a compound that, directly or indirectly, increases the prevalence of the mRNA in 
the cell, the ratio of green to red fluorescence will increase. When the drug decreases the 
mRNA prevalence, the ratio will decrease. 

The use of a two-color fluorescence labeling and detection scheme to define 
alterations in gene expression has been described, e.g., in Shena et al. Science 1995, 
270:467-470. An advantage of using cDNA labeled with two different fluorophores is that a 
direct and internally controlled comparison of the mRNA levels corresponding to each 
arrayed gene in two cell states can be made, and variations due to minor differences in 
experimental conditions (e.g., hybridization conditions) will not affect subsequent analyses. 
However, it will be recognized that it is also possible to use cDNA from a single cell, and 
compare, for example, the absolute amount of a particular mRNA in, e.g., a treated and 
untreated cell. 

Preparation ofmicroarrays. Nucleic acid microarrays are known in the art and 
preferably comprise a surface to which probes tnat correspoiiu m ^^u^^^^. 5— 1— — 
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(e.g., cDNAs, mRNAs, cRNAs, polypeptides, and fragments thereof), can be specifically 
hybridized or bound at a known position. In one embodiment, the microarray is an array in 
which each position represents a discrete binding site for a product encoded by a gene (e.g. , a 
protein or RNA), and in which binding sites are present for products of most or ahnost all of 
5 the genes in the organism's genome. M a preferred embodiment, the "binding site" 

(hereinafter, "site") is a nucleic acid or nucleic acid analogue to which a particular cognate 
cDNA or cRNA can specifically hybridize. The nucleic acid or analogue of the binding site 
can be, e.g., a synthetic oligomer, a full-length cDNA, a less-than fiiU-length cDNA, or a 
gene fragment. 

10 Although in a preferred embodiment the microarray contains binding sites for 

products of all or ahnost all genes in the target organism's genome, such comprehensiveness 
is not necessarily required. Usually the microarray will have binding sites corresponding to 
at least about 50% of the genes in the genome, often at least about 75%, more often at least 
about 85%, even more often more than about 90%, and most often at least about 99%. 
1 5 Preferably, the microarray has binding sites for genes relevant to the action of a drug of 
interest. A "gene" is identified as a segment of DNA containing an open reading frame 
(ORF) of preferably at least 50, 75, or 99 amino acids from which a messenger RNA is 
transcribed in the organism (e.g., if a single cell) or in some cell in a multicellular organism. 
The number of genes in a genome can be estimated from the number of mRNAs expressed 
20 by the organism, or by extrapolation from a well-characterized portion of the genome. When 
the genome of the organism of interest has been sequenced, the number of ORFs can be 
detennined and mRNA coding regions identified by analysis of the DNA sequence. 

Preparing nucleic acids for microarrays. As noted above, the "binding site" to which 
25 a particular cognate cDNA specifically hybridizes is usually a nucleic acid or nucleic acid 
analogue attached at that binding site. In one embodiment, the binding sites of the microarray 
are DNA polynucleotides corresponding to at least a portion of each gene in an organism's 
genome. These DNAs can be obtained by, e.g , polymerase chain reaction (PCR) 
amplification of gene segments from genomic DNA, cDNA {e.g., by RT-PCR), or cloned 
30 sequences. PCR primers are chosen, based on the known sequence of the genes or cDNA, 
that result in amphfication of unique fragments (i.e. fragments that do not share more than 10 
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bases of contiguous identical sequence with any other fragment on the microarray). 
Computer programs are useful in the design of primers with the required specificity and 
optimal amphfication properties. See, e.g., Oligo version 5.0 (National Biosciences). In the 
case of binding sites corresponding to very long genes, it will sometimes be desirable to 
amplify segments near the 3' end of the gene so that when oligo-dT primed cDNA probes are 
hybridized to the microarray, less-than-fiiU length probes will bind efficiently. Typically 
each gene fragment on the microarray will be between about 50 bp and about 2000 bp, more 
typically between about 100 bp and about 1000 bp, and usually between about 300 bp and 
about 800 bp in length. PGR methods are well known and are described, for example, in 
Innis et al, eds., 1990, PCR Protocols: A Guide to Methods and Applications, Academic 
Press hic. San Diego, CA. It will be apparent that computer controlled robotic systems are 
usefiil for isolating and amplifying nucleic acids. 

An alternative means for generating the nucleic acid for the microarray is by 
synthesis of synthetic polynucleotides or oUgonucleotides, e.g., using N-phosphonate or 
phosphoramidite chemistries (Froehler et al. Nucleic Acid Res. 1986, 14:5399-5407; 
McBride et al.. Tetrahedron Lett. 1983, 24:245-248). Synthetic sequences are between about 
15 and about 500 bases in length, more typically between about 20 and about 50 bases. In 
some embodiments, synthetic nucleic acids include non-natural bases, e.g., inosine. As noted 
above, nucleic acid analogues may be used as binding sites for hybridization. An example of 
a suitable nucleic acid analogue is peptide nucleic acid (see, for example, Egholm et al.. 
Nature 1993, 365:566-568. See, also, U.S. Patent No. 5,539,083). 

In an alternative embodiment, the binding (hybridization) sites are made from 
plasmid or phage clones of genes, cDNAs (e.g., expressed sequence tags), or inserts 
therefrom (Nguyen et al.. Genomics 1995, 29:207-209). In yet another embodiment, tiie 
polynucleotide of the binding sites is RNA. 

Attaching nucleic acids to the solid surface. The nucleic acids or analogues are 
attached to a soHd support, which may be made from glass, plastic (e.g., polypropylene, 
nylon), polyacrylamide, nifrocellulose, or other materials. A preferred method for attaching 
the nucleic acids to a surface is by printing on glass plates, as is described generally by 
Schena et al.. Science 1995, 270:467-470. This method is especially usefixl for preparing 

-35- 

{M:\3235\100m087us2\00058277.DOC IMfflBmilllllllllilUlllillllll } 



microarrays of cDNA. See also DeRisi et al. Nature Genetics 1996, 14:457-460; Shalon et 
al. Genome Res. 1996. 6:639-645; and Schena et al, Proc. Natl Acad. ScL USA 1995. 
93:10539-11286. 

A second preferred method for making microarrays is by making high-density 
5 oligonucleotide arrays. Techniques are known for producing arrays containing thousands of 
oligonucleotides complementary to defined sequences, at defined locations on a surface usmg 
photolithographic techniques for synthesis in situ (see, Fodor et al, Science 1991, 251:767- 
773- Pease et al, Proc. Natl Acad. Sci. USA 1994. 91:5022-5026; Lockhart et al. Nature 
Biotech 1996,14:1675. See. also, U.S. Patent Nos. 5.578,832; 5,556,752; and 5,510,270) 
10 or other methods for rapid synthesis and deposition of defined oligonucleotides (Blanchard et 
al. Biosensors & Bioelectronics 1996, 11:687-90). When these methods are used, 
oligonucleotides (e.g., 20-mers) of known sequence are synthesized directly on a surface 
such as a derivatized glass slide. Usually, the array produced is redundant, with several 
oligonucleotide molecules per RNA. OUgonucleotide probes can be chosen to detect 

1 5 alternatively spliced mRNAs. 

Other methods for making microarrays. e.g.. by masking (Maskos and Southern, Nuc. 
AcidsRes. 1992, 20:1679-1684). may also be used. In principal, any type of array, for 
example, dot blots on a nylon hybridization membrane (see, Sambrook et al. Molecular 
Cloning-A Laboratory Manual (2nd Ed.), Vol. 1-3, Cold Spring Harbor Laboratory, Cold 

20 Spring Harbor, N.Y., 1989), could be used, although, as will be recognized by those of skill 
m the art. very small arrays will be preferred because hybridization volumes will be smaller. 

Generating labeled probes. Methods for preparing total and poly(A)^ RNA are well 
known and are described generally in Sambrook et al, supra. In one embodiment, RNA is 

25 extracted from cells of the various types of interest in this invention using guanidinium 
thiocyanate lysis followed by CsCl centrifiigation (Chirgwin et al , Biochemistry 1979, 
18:5294-5299). Poly(A)^ RNA is selected by selection with oligo-dT cellulose (see 
Sambrook et al, supra). Cells of interest may include, but are not limited to, wild-type cells, 
drug-exposed wild-type cells, modified cells, and drug-exposed modified cells. 

30 Labeled cDNA is prepared from mRNA by oUgo dT-primed or random-primed 

reverse transcription, both of which are well known in the art (see. for example. Klug & 
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Berger Methods Enzymol 1987,152:316-325). Reverse transcription may be carried out m 
the presence of a dNTP conjugated to a detectable label, most preferably a fluorescently 
labeled dNTP. Alternatively, isolated mRNA can be converted to labeled antisense RNA 
synthesized by in vitro transcription of double-stranded cDNA in the presence of labeled 

5 NTPs (Lockhart et al. Nature Biotech. 1996, 14:1675). In alternative embodiments, the 
CDNA or RNA probe can be synthesized in the absence of detectable label and may be 
labeled subsequently, e.g. , by incorporating biotinylated dNTPs or NTP, or some similar 
means {e.g., photo-cross-linking a psoralen derivative of biotin to RNAs), followed by 
addition of labeled streptavidin {e.g., phycoerythrin-conjugated streptavidin) or the 

10 equivalent. 

When fluorescently-labeled probes are used, many suitable fluorophores are known, 
including fluorescein, Ussamine, phycoerythrin, rhodamine (Perkin Elmer Cetus), Cy2, Cy3, 
Cy3 5 Cy5 Cy5.5, Cy7, FluorX (Amersham) and others (see, e.g., Kricka, 1992, 
Nonilotopic DNA Probe Techniques, Academic Press San Diego, CA). It will be appreciated 
15 that pairs of fluorophores are chosen that have distinct emission spectra so that they can be 
easily distinguished. 

In another embodiment, a label other than a fluorescent label is used. For example, a 
radioactive label, or a pair of radioactive labels with distinct emission spectra, can be used 
(see Zhao era/., Gene 1995, 156:207; Pietuer a/., Ge«omei?e.. 1996,6:492). However, 
20 because of scattering of radioactive particles, and the consequent requirement for widely 
spaced binding sites, use of radioisotopes is a less-preferred embodiment. 

hi one embodiment, labeled cDNA is synthesized by incubating a mixture containing 
0.5 mM dGTP, dATP and dCTP plus 0.1 mM dTTP plus fluorescent deoxyribonucleotides 
(e g , 0.1 mM Rhodamine 1 10 UTP (Perken Elmer Cetus) or 0.1 mM Cy3 dUTP 
25 (Am'ersham)) with reverse transcriptase (e.g.. SuperScriptTM. H, LTI Inc.) at 42 °C. for 60 



min. 

Hybridization to microarrays. Nucleic acid hybridization and wash conditions are 
chosen so that the probe "specifically binds" or "specifically hybridizes" to a specific array 
30 site, i.e. , the probe hybridizes, duplexes or binds to a sequence array site with a 
complementary nucleic acid sequence but does not hybridize to a site with a non- 
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complementary nucleic acid sequence. As used herein, one polynucleotide sequence is 
considered complementary to another when, if the shorter of the polynucleotides is less than 
or equal to 25 bases, there are no mismatches using standard base-pairing rules or, if the 
shorter of the polynucleotides is longer than 25 bases, there is no more than a 5% mismatch. 
Preferably, the polynucleotides are perfectly complementary (no mismatches). It can easily 
be demonstrated that specific hybridization conditions result in specific hybridization by 
carrying out a hybridization assay including negative controls (see, e.g., Shalon et al, supra; 

and Chee et al, supra). 

Optimal hybridization conditions will depend on the length (e.g., oligomer versus 
polynucleotide greater than 200 bases) and type {e.g., RNA, DNA, PNA) of labeled probe 
and immobilized polynucleotide or oUgonucleotide. General parameters for specific {i.e., 
stringent) hybridization conditions for nucleic acids are described in the definitions provided 
in Section 5.1 , supra. When cDNA microarrays, such as those described by Schena et al. are 
used, typical hybridization conditions are hybridization in 5x SSC plus 0.2% SDS at 65 °C 
for 4 hours, followed by washes at 25 °C in low stringency wash buffer {e.g., Ix SSC plus 
0.2% SDS) followed by 10 minutes at 25 °C in high stringency wash buffer (O.lx SSC plus 
0.2% SDS). See, Shena et al, Proc. Natl Acad. Sci. USA 1996, 93:10614). Useful 
hybridization conditions are also provided in, e.g., Tijessen, 1993, Hybridization With 
Nucleic Acid Probes, Elsevier Science Publishers B.V. See, also, Kricka, 1992, Nonisotopic 
DNA Probe Techniques, Academic Press, San Diego, CA. 

Signal detection and analysis. When fiuorescently labeled probes are used, the 
fiuorescence emissions at each site of a transcript array can be preferably detected by 
scanning confocal laser microscopy, hi one embodiment, a separate scan, using the 
appropriate excitation line, is carried out for each of the two fluorophores used. Alternatively, 
a laser can be used that allows simultaneous specimen illumination at wavelengths specific to 
the two fluorophores and emissions fi-om the two fluorophores can be analyzed 
simultaneously (see. Shalon et al. Genome Research 1996, 6:639-645). In a preferred 
embodiment, the arrays are scanned with a laser fluorescent scanner with a computer 
controlled X-Y stage and a microscope objective. Sequential excitation of the two 
fluorophores is achieved with a multi-line, mixed gas laser and the emitted light is split by 

-38- 

{M:\3235\100m087us2\00058277.DOC I 



wavetogtt. and detected with two photomulUpUer tubes. Fluorescence laser scannmg devices 
are described in Schena « Genome Res. 1996. 6:639-645 and in other references Ced 
herein. Alternatively, the fiber-optic bundle describ«i by Ferguson <./., Namre Biotech. 
1996, 14:1681-1684, may be used to monitor mRNA abundance levels at a large number of 

sites simultaneously. 

Signals are recorded and, in a preferred embodiment, analyzed by computer, e.g., 
using a 12 bit analog to digital board. In one embodiment the scamied image is despeckled 
using a graphics program (e.g., Hijaak Graphics Suite) and then analyzed using an image 
gridding program that creates a spreadsheet of the average hybridization at each wavelength 
at each site. If necessary, an experimentally determined correction for "cross talk" (or 
overlap) between the channels for the two fluors may be made. For any particular 
hybridization site on the transcript array, a ratio of the emission of the two fluorophores can 
be calculated. The ratio is independent of the absolute expression level of the cognate gene, 
but is useful for genes whose expression is significantly modulated, e.g., by administenng a 
drug, drug-candidate or other compound, or by any other tested event. 

' m one preferred embodiment of the invention, the relative abundance of an mRNA m 
two cells or ceU lines tested (e.g., in a treated verses untreated cell) may be scored as 
perturi,ed {i.e., where the abundance is different in the two sources of mRNA tested) or as 
not perturbed (i.e., where the relative abundance in the two sources is the same or is 
unchanged). Preferably, the difference is scored as perturbed if the difference between the 
two sources of RNA of at least a factor of about 25% (i.e., RNA from one sources is about 
25% more abundant than in the other source), more preferably about 50%. Still more 
preferably, the RNA may be scored as perturbed when the difference between the two 
sources of RNA is at least about a factor of two. Indeed, the difference in abundance 
between the two sources may be by a factor of three, of five, or more. 

m other embodiments, it may be advantageous to also determine the magnitude of the 
perturbation. This may be done, as noted above, by calculating the ratio of the emission of 
the two fluorophores used for differential labeling, or by analogous methods that will be 
readily apparent to those of skill in the art. 
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6.5. Rimnformat ir.s and Statistics 

Those skilled in the art will readily appreciate that the MPHTS assays of this 
invention will, at least in preferred embodiments, track a large amount of data from many 
sources including, e.g., expression levels for a large number of different genes in a variety of 
5 different cell and tissue types and under a variety of different conditions. The invention 
therefore preferably makes use of methods in bioinformatics and statistical analysis to 
integrate such data. Such analysis tools include, for example, clustering and class 
partitioning algorithms that enable a user to summarize and visualize effects of multiple 
variables on relationships within a data set. In a particularly preferred embodiment, the 
10 MPHTS methods of this invention make use of a statistical analysis tool referred to as 

'Trincipal Component Analysis" or "PCA". The technique is well known in the art and may 
be implemented, e.g., using commercially available software such as the Partek suite of 
pattern recognition tools (Partek Inc., St. Charles, Minnesota). 

By PCA analysis of gene expression data from different brain areas and disease 
15 states, a user is able to readily identify whether the major source or sources of variance 
within the data set are correlated with the particular cells or tissue and/or whether such 
variance is correlated with a neuropsychiafric disorder of interest. An exemplary figure 
depicting this analysis is set forth here, in FIG. 2. Those skilled in the art will readily 
appreciate and/or be able to select appropriate cutoffs {e.g., a maximum significant p-value) 

20 for use in such methods. 

Statistically significant changes in gene expression may also be identified by 
coordinately regulated genes in distinct pathways, as well as coordinate changes of multiple 
genes within a common pathway {e.g., genes involved in a common metaboHc pathway or 
process). These provide an aggregate level of statistical significance that far exceeds the 
25 statistical significance obtained for the genes individually. 

hi preferred embodiments, RNA extraction and/or hybridization experiments are 
repeated at least once, and more preferably multiple times for each sample to assure 
statistically robust and reproducible results. Changes in gene expression that appear to be 
statistically significant may also be confirmed by an independent experimental technique 
30 such as real-time polymerase chain reaction (RT-PCR), quantitative in ./Yuhybridization, 
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immunohistochemistry and functional assays of the translated protein(s), all of which are 
well known and routinely used in the art. 

6.6. T Tses of Gene <iipnatiires in MPHTS 

5 Once genes signatures or "efficacy genes" for a particular disorder have been 

identified and/or selected, they may be readily used, e.g., in screening assays to identify 
promising therapeutic compounds that may be useful for diagnosing and/or treating such a 
disorder. In particular, a candidate therapeutic compound may be identified in screening 
assays of the invention by identifying compounds that produce changes in the expression of 

10 signature genes that are similar to the changes observed in the ECS model, demonstrated in 
the exmples, infra. Such changes may be identified qualitatively, but are more preferably 
identified quantitatively; for example, by assigning a MPHTS "value" for each compound 

tested in the screening assay. 

As an example, and not by way of limitation, such an MPHTS value may simply be 
15 the sum of changes in each signature gene's expression observed for a test compound in the 
screening assay. Preferably, these changes in the signature genes' expression levels are 
normalized as apercentage of the "optimal" change in each gene's expression. As used here, 
the change in expression of an signature gene is said to be "optimal" when it is 
approximately equal to the change in expression associated with a therapeutic benefit. 
20 Optionally, the change in each signature gene's expression can also be weighted, e.g., by the 
signature gene's score (as determined, e.g., using the algorithm set forth in Section 8.2, 
infra)- calculation of such a value can be easily represented mathematically by the 
equation: 

V = Yu^iEi (Equation!) 

i 

25 In Equation 1, Fis the MPHTS "score" calculated for a test compound in an MPHTS 

assay is the measured change in the expression of gene i in cells contacted with the test 
compound compared to the expression in cells that are not contacted with the test compound. 
As noted above, 5, will preferably be normalized to the "optimal" change associated with a 
desired therapeutic effect. For example, can be expressed as the percentage or firaction of 

30 optimal change, w, indicates the score tor me signature gcuc y......^^ ^^^^^ . 
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is obtained or derived from the score value calculated for gene / (e.g., according to the 
algorithm set forth in Section 8.2) and is converted to a percentage of the average score value 
for the signature genes that comprise the entire set used for drug screening. 



7. EXAMPLES 

The present invention is also described and demonstrated by way of the foUowmg 
examples. However, the use of these and other examples anywhere in the specification is 
illustrative only and in no way limits the scope and meaning of the invention or of any 
exempUfied term. Likewise, the invention is not limited to any particular preferred 
embodiments described here. Indeed, many modifications and variations of the invention 
may be apparent to those skilled in the art upon reading this specification, and such vanations 
can be made without departing the invention in spirit or in scope. The invention is therefore 
to be limited only by the terms of the appended claims along with the full scope of 
equivalents to which those claims are entitled. 

7 1 Gex^e SienaturftR nf Electros hnck Therapy 

a« Modeled b v indent Electr "''^"vi'1sive Seizures 

7.1.1 Materials and Metiiods 
Subjects. Male, Sprague-Dawley rats (ACE Animals; Boyertwon, PA) were 2-6 
months of age and weighted 200-350 g at the start of each experiment. Animals were housed 
2 per cage, understandard conditions, including ad libitum access to food and water m a 
12:12 hour Ught/dark vivarium. 

ECS Delivery. Animals were handled for 3 days prior to the begimiing of the 
experiment for them to become adjusted to the seizure induction area and to apphcation of 
the ear-clip electrodes. The seizure induction or sham-seizure induction areas consisted of 
plastic cages with fresh bedding material. Each cage was placed on a lab bench and a pair of 
electrodes was hung over the top of each cage. One set of electrodes was comiected to the 

ECT Unit (Ugo Basile # 57800, Comeno, itaiyj lor use wiui uxw ^^i> ^ 

set (for control animals) was discomiected from the ECT unit such that no current could be 
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delivered. Cage mates were run concurrently: one animal was randomly assigned to the 
shock group and the other to the control group. Both animals were removed from their home 
cage and placed in the appropriate testing cage. Each set of ear-clip electrodes were 
dampened with saline and attached to the deep region of the outer ear. The shock, a square 
wave pulse (pulse width = 0.7 ms, 50 Hz) of 50 mA (Madsen et al.,Biol. Psychiatry 2000, 
47:1043-1049) was administered for a duration for 1.0 second to provide a tonic phase 
characterized by extension of all limbs and forward head extension, lasting for 10- 
15 seconds. Animals, both ECS and control, were returned to their home cage 10 minutes 
after seizure induction. Acute ECS animals (N = 12) were given a single shock. Chronic 
ECS animals (N = 12) were given one ECS per day for 10 consecutive days. Control group 
rats (N = 12 acute and N = 12 chronic) followed this protocol. 

Visual observation of behaviors, side effects. Visual observations were sufficient to 
confirm the presence of clonic-tonic seizures immediately following the shock. Side effects 
during the post-ictal phase, which lasted about 10 minutes, were best characterized by 
heightened responsivity to auditory and tactile stimulation and decreased motor activity. 
General behavioral observations made in the home cage for up to 6 hours post-shock 
administration revealed sedation but no signs of recurrent seizures. 

Brain area dissections. Animals were sacrificed by decapitation without anesthesia 
approximately 3-5 hours after the final ECS treatment. The brain was removed from the 
calvarium and immersed in crushed ice for 3-5 minutes. Coronal slabs of 2 mm thickness 
were prepared, from which the frontal cortex and hippocampus were dissected. The brain 
regions were weighed, frozen on dry ice, and stored at -80 "C. 

RNA Extraction. Total RNA was extracted from about 70-100 mg of frozen rat tissue 
with 0.7-1.0 ml of TRIzol reagent (hivitrogen; Carlsbad, CA) as per the manufacturer's 
instructions. The tissue was homogenized with a Polytron (Fisher) for 30 seconds. The 
homogenate was cleared by a 10 minute centrifiigation at 10,000 x g, and 0.2 ml of 
chloroform was added per ml of TRIzol. The samples were vigorously shaken for 
20 seconds and incubated on ice for 10 minutes. The aqueous phase was separated by 
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centrifugation at 10,000 x g for 10 rr^inutes, allowed to precipitate for 10 minutes, and the 
preapitatecollectedbycentrifugationatl2,000xgforl0minutes.Tl.e^^^^^^ 

with 700/0 ethanol collected by a brief centrifugation, air dried, and resuspended m HA The 
RNAwas further purified using an RNeasycolunmCQiagen; Valencia, CA). Thepunfied 
RNA was quantified by a UV spectrophotometer, and RNA quality determined by capUlary 
electrophoresis on an Agilent Bioanalyzer. 

Gene expression profiling. RNA (10 ^g) was converted to double-stranded cDNA 
following priming with an oligo-dT-T7 primer (Operon). The resultant cDNA was punfied 
by phenol-chloroform extraction, ethanol-precipitated, and resuspended in H,0. The punfied 
cDNA was subjected to in vitro transcription using T7 polymerase in the presence of 
biotinylated UTP and CTP (Enzo Ufe Sciences, Farmingdale, NY). The resultant aRNA was 
purified with an RNeasy column (Qiagen, Valencia, CA), eluted in HA and quantified usmg 
a UV spectrophotometer. aRNA (15 ,g) was fragmented following the Affymetrix protocol. 
> The quality of RNA fragmentation was checked on the Agilent Bioanalyzer, and 30 ^1 was 
added to 270 ,il of hybridization buffer and hybridized to the Affymetrix U34A Rat Genome 
GeneCWp®. After 16 hours of hybridization at 45 °C the GeneChip was washed, stained, 
and scamied according to the standard Affymetrix protocol. 

0 StatisUcal analysis. Each gene on Affymetrix U34A gene chips is represented by one 

or more of about 8800 "probe sets" - i.e. , sets of oligonucleotides that are derived from the 
mRNA sequence of that particular gene. Following hybridization and scam^ing, each probe 
set was assigned a '^present" or "absent" call and an overall fluorescent signal intensity as 
determined by the Affymefrix Microarray Suite Software 5.0 following global normahzatxon. 

.5 A filtering of all probe sets was conducted to remove those genes that were clearly not 

expressed in the tissue, and such probe sets were not analyzed. Specifically, a probe set for a 
particular gene was retained for all samples only if it contained a "present" call for at least 
one sample, and only if the mean fluorescence intensity for either the control or treatment 
group was greater than 100. This procedure filtered out about 50% of all probe sets, 

30 consistentwiththeknownexpressionofapproximately40-50o/oofratgeneswithinbram 

tissue. Absent or present calls are assigned by Affymetrix software on the basis of a 
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comparison between the intensities of perfect oUgonucleotide matches and smgle 
oligonucleotide mismatches. THe finding that mismatch intensities are not sigmficantly 

lower than those for the perfect match oligonucleotides does not necessarily imply that the 
gene is not expressed, as it may be a result from cross-hybridization of other mRNA species, 
5 other technical artifacts, or measurement noise. Therefore, present calls were counted m a 
liberal mamier to avoid discarding probe sets where there was some indication that the gene 
may in fact be expressed. The alteration of many genes in more than one tissue or after both 
durations of ECS, and the confirmation of many gene changes by RT-PCR, confirms that 
measurements of some genes with relatively few present calls would be detected m 
10 independent tissues and by an independent method. 

A principal components analysis (PCA) was performed using Partek Pro 2000 
software package (Partek, Inc.). The analysis treated arrays as "objects" and probe sets as 
"variables" to detect outlier arrays. PCA plots mapped the distidbution of arrays wrthm a 3- 
dimensional space defined by the three greatest contributors to the variance in gene 
15 expressionvalues. Tl^epositionof an array on such plot refiects the overall expression 
profile of the sample as measured by the array. The arrays that appeared as outliers on the 
plots were removed from fiirther consideration. Such PCA outliers are usually associated 
with a technical problem, such as marginally acceptable RNA quality, regional blotches m 
hybridization signals, or weak overall hybridization on the corresponding gene chip. For 
20 each of the four experiments, (acute and chronic ECS, frontal cortex and hippocampus), at 
least seven, but typically nine or all 10, of the rat samples from each conti-ol or tireatinent 
group were retained for fiirther analysis. 

The ability of acute or chronic ECS to alter gene expression versus that of sham- 
treated controls was determined for each tissue type. The magnitude of expression level 
25 change was estimated as the ratio of the means between the tireated and confrol groups. 
Probe sets were considered as "hits" only when the average normalized signal intensity 
exceeded 100 in either of the groups, the ratio of means exceeded 1.5 or were less than 0.67 
and the p value of the two-tailed Welch t test was less than 0.05. 

In a multiple testing setting, the p value of 0.05 does not necessarily mean that the 
30 probability of type 1 error is really 0.05, unless the Bonferroni's correction has been applied. 
Bonferroni's correction was not applied in these experiments, however, because it would 
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,esul.mmany.ype2«Tors andteteofmos. valid hits. The uncorrecedpvalueofa 
statistical .^tisftereforeu^fulhe^primarilyasameasure of diffe^ncebetweendataseu. 

According to a «cen. study (Xu & Li. BioinformaUcs 2m. 19:1284-1289). a parametric p 
value „,ay outperfonn a pennutational one as a distance metric. Additionally, the 1.5 rauo of 
5 means cut-offused here increased the verification rate of gene expression changes. 

GeneSpring® analysis software (SiUcon Genetics) was used to evaluate the overlap 
between the gene changes identified m the four comparisons. Since about 3,500-4,400 
genes were pre«>t in each sample, a significant number of false positives was expected. To 
help remove false positives, and also to increase confidence in selecting gene tats for RT- 
10 PCR confirmatton, each statistically significant gene change was ranked objectively by a 
biostatistical ECS gene selection algorithm (set forth in the appendix at Section 8.2. .nfra) 
according to its statistical and biological significance, and by the robustness of experimental 
change, including whether it was altered in one or both tissues and after acute, chrome, or 
both ECS durations. 



15 



Real-time polymerase chain reaction assay (RT-PCR). Total RNA (2 ^g) were 
subjected to DNAse treatmer^t in a 10 ,1 reaction containing 1 ^ lOX DNAse I reaction 
buffer, and 1 .IDNAse I (Invitrogen, Carlsbad, CA). The reaction was carried out at room 
temperature for 10 minutes. OnettlofEDTA(25mM)and miofoligo(dT),2-,8mer(0.5 
20 ^g/ml, Invitrogen, Carlsbad, CA) were added to DNAse reaction and heated to 70 °C for 1 5 
minutes in a water bath to simultaneously inactivate the DNAse I activity and elimmate RNA 
secondary structure to allow oligo dT-poly A annealing. The sample was placed on ice for 2 
minutes and collected by brief centrifugation. The RNA in the sample was rexerse- 
transcribed into cDNA by the addition of 8 ^1 of master mix containing 4 ^1 of 5X first 
25 strand buffer, 2 ^1 DTT (0.1 M), 1 ^1 dNTP's (10 mM each), and 1 ^1 Superscript II 
(200U/^l) (hivitrogen, Carlsbad, CA), and incubated at 42 °C for 45 minutes. The RT 
reaction was diluted 10 fold with dHjO and stored at 4 °C. 

For each specific gene tested, diluted cDNA (5 ^l) was added to 45 ^ PCR reaction 
mixture that contained 25 ^tl of 2X SYBR^ Green PCR Master Mix (Applied Biosystems. 
30 Warrington, UK) and 25 picomoie ot eacn lorwaru aim icvo... px....... ^ --r 
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subjected to 40 cycles of Real Time PCR (ABI 7900, Applied Biosystems, Wamngton UK) 
where fluorescence is measured several times during each cycle of 2-^ PCR alternating 
between about 95 T for 15 seconds and 60 T for 1 minute. The thr«hold cycle (Ct), or 
point at which signal fluorescence exceeds background, for each sample for each gene was 
5 eompared^-astandardcurvetodetenninearelativeexpressionvalue. The standard curve 
was generated by real time PCR analysis of five, 1 0-fold serial dilutions of a cDNA 
generated 6om the RNA of rat cortical stem cells. This method allows a relative 
comparison between samples. The expression value of each gene was nomtalized to the 
relative amount of GAPDH expressed in that sample to calculate a relative amount of 
,0 transcriptpresentforeachgene. The nonnalizrf expression values for all control and tteated 
samples were averaged and an average fold change determined. A Smdenfs . test was 
conducted between the normalized relative expression values for each individual control and 
treated samples to determine statistical relevance. 

15 7.1.2 Results . 

Four sets of microarray hybridizations to evaluate the effects of acute or chrome ECS 
treamtentongeneexpressionintheftontalcortexandhippocampus. Plots of the mean 
abundances of each probe set measured with the Aff^etrix U34A rodent chip are illustrated 
in Figures 1 A-ID the frontal cortex (Figure, lA and IB) and hippocampus (Figures IC and 
20 ID) of rats that wUh exposed to either acute (Figure, lA and IQ or chronic (Figure, IB 
and ID) ECS, compared to sham control groups that received no ECS. At gene abundances 
above 100 each scatty plot of mean intensities diplayed minimal skewing about the umty 
Une The di^bution of points was generally quite close to the unity line, and stat,st,cally 
significant outliers (illustrated by white dots in Figure, lA-lD) were frequently weU outside 
25 dte main duster, of genes. Principal components analyri, (PCS) revealeded that the frontal 
cortex and hippocampus transcriptome segregated into distinct duster, similar to that 
reported for different human CNS btair, region, (Pafeyman e, al.. Current Drug Targets 
2002 1 ■205-216). Outlier samples in the PCA analysU could be readily distingmshed, and 
eithe^ one sample, two samples or no sample group were removed fi^m the statistical 
30 analysU so that the final n was 8-10 per group. 
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Data quality evaluation. In order to first determine the number of samples per group 
r^eeded to rdiably identify differentially expressed genes, subgroups of 2 to 10 samples were 
randomly created from the sham- or ECS-treated groups (or fi:om these two groups 
combined) of gene expression data in the acute hippocampus. Welch t tests were performed 
5 for all gene changes comparing subgroups within or between the sham control and ECS 
groups as well as between "mixed" subgroups, and this was repeated for 100 different 
comparisons. Each time, the number of "hits" passing p value cut-off of 0.01 was counted, 
and these data are plotted in Figure 2. Statistical theory predicted that about 40 genes are 
expected to appear as hits due to random technical and biological variations fi-om the 3,500 to 
10 4400detectedgenes,independentofanyrealECStreatmenteffect. Comparisons between 
subgroups of samples drawn from within a group of animals or from the "mixed" ammals 
produced numbers of significant gene changes that were mostly below this theoretical chance 
level and gradually approached it with increases in subgroup size, as can be seen from visual 
inspection of Figure 2. In contrast, the number of significantly changed genes in the analysis 
15 of sham- versus ECS-treated animals increased steeply from chance level at group size of 2. 
to 15-fold above chance at group sizes of 10. Tliis analysis demonstrates that there is 
actually a highly significant effect of the ECS shock treatment on gene expression measured 
by the microarrays. The analysis also demonstrates that, while greater numbers of ammal per 
group increase statistical power, only 8-10 animals per group are necessary to reliably detect 
20 large numbers of changes in gene expression levels. 

Identification and ranking of differentially expressed genes. A total of 135 unique 
genes represented by more than 200 probe sets, were found to be differentially expressed m 
ECS-ireated rats compared to controls in at least one of the four different ECS experiments. 

25 These genes are Hsted in the Appendix, infra (see. in particular, at Section 8.3) along with the 
Accession number for their nucleotide sequence on the GenBank database. THe Table m 
Section 8 3 also provides the expression ratio and p-value determined for each of these genes 
in the microarray analysis, as well as the "score" determined using the algorithm set forth at 
Table X, supra. Many of these genes, whose names are shaded in the table at Section 8.2, 

30 have not been previously associated with either ECS or ECT. 
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A Venn diagram representing the distribution of changes among the 135 genes in 
different tissues and ECS treatment groups is shown in Figure 3. Forty-four genes were 
identified that had significant differential expression in multiple tissues and/or m mulUple 
ECS treatment groups. The most dramatic response to ECS treatment occurred m the 
hippocampus, where 87 genes were detected as differentially regulated following a smgle, 
acute shock. Chronic ECS treatment resulted in fewer gene changes in the hippocampus (64 
compared to 87), about half of which were distinct fi:om those that changes after a smgle 
ECS The frontal cortex was much less responsive than the hippocampus, with 36 genes 
identified as differentially regulated following acute treatment, and 36 genes identified as 
differentially regulated by chronic ECS treatment. Twelve genes that were regulated m the 
same direction in the hippocampus by ECS and following exercise (Moltem et al., Eur. J. of 
Neurosci. 2002, 16:107-1 16; Tong et al.,Neurobio. of Disease 2001, 8:1046-1056) are 
indicated in Section 8.3 by an asterix (*) next to the gene name. A supplemental list of genes 
with p values less than 0.005 and having differential expressing changes less than 1.5 fold 
was also identified. THese genes are listed separately, infra, in Section 8.4 infra. As with the 
genes listed in Section 8.3, many of the genes listed in Section 8.4 have not been previously 
associated with either ECS or ECT.niese genes' names are shaded in the tables at Sections 



20 



25 



8.3-8.4. 



Many genes are represented by more than one probe set on the Affymetrix RG-34A 
chip Among genes exhibiting significant change due to ECS treatment, 13 were identified 
by multiple probe sets. These genes are indicated in Tables of Sections 8.3 and 8.4 by 
italicized type in those tables. Five of these thirteen genes were selected to determine if all of 
the probe sets for each gene were in agreement. The numbers of probe sets reportmg 
significant and co-directional change due to ECS treatments were as follows: five out of 
seven for BDNF; two out of two for Cox-2; two out of four for Jmi, two out of three for 
TIEG and 1-2 out of four (depending on the treatment and tissue) for catalase. The 
remaining five probe sets for BDNF, Jun and TIEG, as well as one of the catalase probe sets, 
produced only absent calls and/or less than 100 abundance units in all rat brain samples. 
These findings suggests that some of these probe sets simply failed to measure their target 
30 mRNA, and that ECS actually did affect at least four out of these five genes. 
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Twenty-six of the identified genes in Section 8.3 were found to be similarly up- or 
down-regulated in both the hippocampus and frontal cortex (Figure 3), indicating that many 
brain-wide changes following shock therapy may occur independent of cellular heterogeneity 
of these different regions. M addition, 39 gene were found to change after both acute and 
chronic ECS treatments. Of the genes that were either up- or down-regulated in more than 
one tissue type and/or after both acute and chronic ECS therapies, one gene changed in 
different directions. This gene, solute carrier number 3, member 1, decreased expression in 
both the frontal cortex and hippocampus following chronic ECS, but increased expression in 
the hippocampus after acute ECS. The expression of five genes in particular was similarly 
affected by ECS in every experiment (labeled as "all" in the Table of Section 8.3) - i.e., 
following both chronic and acute ECS treatment in both the hippocampus and frontal cortex. 
These five genes are: brain derived neurotrophic factor (BDNF), prostaglandin- 
endoperoxide synthase 2 (COX-2), neuronal activity-regulated pentraxin (Narp), TGF^- 
inducible early growth response, and tissue inhibitor of metalloproteinase 1 (TIMP-1). Each 
of these genes or their biological pathways have been previously reported to change 
individually following either ECS treatment or seizure in rats (see, in particular, Zetterstrom 
etai, Brain Res. Mol. Brain Res. 1998,57:106-110; mhuya et al., J. Neurosci. 1995, 
15:7539-7547; mdKzshimoto, Brain Res. 1998,804:212-223). However, the combined 
changes in expression of these genes has not heretofore been demonstrated. 

About two-thirds (i.e., 92 out of 135) of genes listed in Section 8.3 are changed their 
expression in only a single experiment. Hence, acute and chronic ECS can influence gene 
expression in different ways, and genes in the hippocampus and frontal cortex often respond 
differently to ECS treatment. 

Verification of gene signatures by RT-PCR. The measurement of expression levels 
for a large number of genes can result in a significant number of false positive results 
(Zhang, Proc. Nat. Acad. Sci. U.S.A. 2002, 99:12509-1251 1; Zhou et ah, Proa Nat. Acad. 
Sci. U.S.A. 2000, 99:12783-12788; and Xu & Li, Bioinformatics 2003, 19:1284-1289). To 
confirm the validity of gene expression changes in these experiments, the expression levels 
of sixty genes were retested using quantitative real time polymerase chain reaction (RT-PCR) 
using an mdependent sampling of the RNA samples in the microarray analysis. 
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Approximately 1,200 RT-PCR assays were conducted with hippocampal or frontal cortex 
tissues after actue or chronic ECS, focussing on those genes with the largest fold change and 
the best statistical signficance in the microarray analysis. Genes were also selected for 
confirmation by RT-PCR to provide a representation of various protein classes and for their 
biological relevance to depression or its treatment. Hie results of this analysis are presented 
infra, in Section 8.5 of the appendix. Genes were identified with a liberal p < 0.05 cntena 
but required at least a 1 .5-fold increase or decrease to be considered significant. The 
abbreviations "HCA" and "HCC" refer to genes identified in the hippocampus following 
either acute or chronic ECS treatment, respectively. The abbreviation FCA indicates that a 
gene was identified in the prefontal cortex following acute ECS treatment. 

A "first pass" validation of the microarray results indicated a 50% validation rate. 
Principal failures in these experiments were apparently mostly caused by low technical 
performance of individual primer pairs, as evidence by a failure of the ampHfied signal to 
reach asymptotic intensity levels characteristic of genes that are efficiently amplified. For 
those primer pairs that failed, a second and unique primer pair was selected from the same 
region of the franscript to generate an amplicon the same length as the original primer pair. 
Non-validated genes were retested with the new primer pairs.80o/o of these RT-PCR retested 
genes produced the expected asymptotic signal across their amplification cycles, and 
validated results obtained in the microarray analysis. 

The RT-PCR experiments confirmed statistically significant differential expression 
for 43 out of the 60 genes tested {i.e., 63%). hi many instances, the average magnitude of the 
change in expression for a given gene was only 30% between the two groups. However, 
individuals within those different two groups would often differ by a factor of two or more. 
Hence, the RT-PCR procedure was used to quantify differences between individuals that 
exceed differences between the different experimental groups. 

RT-PCR produces a false negative rate of approximately 30% for any large set of 
genes that has been found to be significantly changed on a microarray platform. Only a few 
genes exhibited high, asymptotic intensity levels during RT-PCR while, at the same time, 
failing to replicate microarray results because of low statistic significance. However, these 
nevertheless displayed a magnitude and direction of gene abundance change that were similar 
to those seen with the microarrays. Therefore, it is more likely that greater than 63%. and 
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possible closer to 8O0/0 of the genes detectedby microarray analysis of ECS tissues (i.e., m 
Sections 8.3 and 8.4) represent true positives based on RT-PCR validation. Repeated 
observations of expression changes within different tissues for these genes, confirmation of 
the magnitude and direction of expression changes for 63 to 80«/o of these genes by RT-PCR, 
and the substantial agreement of many observed changes with other literature findmgs all 
confirm that the results reported here are valid. 

7 2 Rplprtion of Gen'^^ fnr KCS Gene Signatures 

To select differentially expressed genes that are most informative and therefore usefiil 
ECS signature, e.g. , for MPHTS type screening assays of the invention, the genes listed 
Sections 8.3 and 8.4, infra, were prioritized for overall significance using the algorithm set 
forth in Section 8.2. This algorithm is based on statistical, biological and experimental 
considerations that account for factors such as a gene's known relationship (if any) with 
bipolar affective disorder and/or unipolar major depression - as well as accounting for the 
magnitude and/or statistical significance (e.g., the p-value) of the gene's differential 
expression in one or more ECS experiments. The algorithm minimizes biases (for example, 
favoring genes that change in one particular region or tissue) that are likely to occur when 
other methods of selection are used. As such, algorithms such as in the one in Section 8.2, 
infra, provide an objective and preferable way to identify informative genes that can be used 

in an ECS gene signature. 

For each gene Usted in Sections 8.3-8.4, a numerical value was assigned for the three 
aspects listed in Section 8.2 - i.e., for statistical considerations (I), for biological 
considerations (H) and for experimental considerations (III). The numerical value(s) for each 
consideration were computed as described in Section 8.2. For biological considerations (II) a 
numerical value was assigned for each gene based on reports fi:om the literature (if 
available). 

The algorithm score thus calculated for each gene is set forth in the far right-hand 
column of the tables at Sections 8.3-8.4. 18 out of the 20 genes most highly ranked by the 
algorithm have been identified by prior investigators to change after ECS, whereas none of 
) the 30 most lowly ranked genes have been previously associated with ECS. Yet, prior 
knowledge about a gene's change in response to ECS accounts for only about 2% of the 
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algorithm's total score, demonstrating that the algorithm is able to independently prioritize 
genes according to their known response in ECS. In contrast, ranking genes according to 
only the level of their expression change and/or its p-value produces a more random 
distribution of the genes that have been previously associated with ECS. 

5 

To better identify differentially expressed genes that are particularly preferable, 
.for use in an ECS gene signature, a biological consideration score (II) was recalculated 
for lach gene - this time using experimental data from homologous genes in human cells 
and/or tissues. Specifically, expression levels of homologous genes were obtained from 
10 "disease signatures" in human tissue samples that had been obtained from individuals 

diagnosed with a neuropsychiatric disorder. In additiona, expression levels of homologous 
genes in "drug signatures" from neuronal cells or animals treated with neuropsychiatric drugs 
{e.g., valproate, clozapine, carbamazepine, and lithium) were also considered. Such "drug 
signatures" and "disease signatures" have been previously described, as have routine 
15 methods and techniques by which they can be obtained. See, in particular, U.S. patent 
application Publication No. 2003/096264 Al . 

Homologos and/or orthologs of the rat genes listed in the appendices, infra, can be 
readily identified, for example by their level of sequence identity to the recited rat nucleic 
acid sequences, or by the level of sequence identity and/or homology to amino acid 
20 sequences that they encode. Alternatively, homologs and orthologs (including those from • 
other species, such as humans and mice) can be identified by hybridization under conditions 
of appropriate stringency, such as those described, supra, in Section 6.1. hi a preferred 
embodiment, appropriate homologs and/or orthologs are identified by routine clustering 
algorithms. For example, the NCBI UniGene database (accessed via the URL 
25 <http://www.ncbi.nlm.nih.gov/UniGene>) groups genes into appropriate clusters of 

homologous sequences from the same and/or different species of organism. See. Schuler, J. 
MoLMed. 1997, 75(10):694-698; Schuler a/., Sc/ence 1996, 274:540-546; andBoyuski& 
Schnler, Nature Genetics 1995, 10:369-371. 

Section 8.6 below provides a list of genes that are differentially expressed in the ECS 
30 experiments (Section 7.1, supra) along with the rat, human and mouse UniGene cluster 
number for each gene and its homologs. The GenBank Accession number for a 
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representative nucleotide sequence (or fragn,ent thereof) for each gene U also provtded 
Overiapping expression data for each gene in human tissue, human cells, and mouse/rat 
tissue is also provided, where available. T.e column "PG. ir gives ,h. algorithm score for 
biolo^ca. considerations (11) calculated using thts additional information. The next column 
("Gronps I. n (Literature) and ffl") provides algorithm scores computed using btologtcal 
eonsideration information ftcm the UteraUtre - U,. the algorithm scor. set fonh m precedmg 
Sections 8.3 and 8.4. Tbe far right-hand column it, Section 8.5 gives the "ECS algonthm 
score" - i.e., the sum of the original score and PGI U. 

The actual number of ECS signature genes used in apatticular appUcatton (..e. for a 
p^cular "gene signature) may vary and generally depends on the need for addtUonal 
information afforded by including more signamre genes, balanced by considerations such as 
the cost and/or efTort involved in addition additional signature genes, e.g., m amtcroa^y or 
PCR assay platfonn. As an example preferred gene signatures of the invention compnse 
nucleic acids cor^sponding to five or more of the genes set forth in Sections 8.3-8.5, 
I„ott.erembodiments,agenesignatureoftheinventionc^comprise„ucle,cac,ds 

corresponding to and/or capable of detecting at least 5. 10. .5. 20, 25, 30. 35, 40, 45. 50, 60, 
70, 80. 90, 100, 1 10, 120, 130. 140. 150 or more of the signature genes set forth m Sections 

8.3-8.5, infra. 

Preferably. signaUire genes that are most infonnative are used. Such genes can be 
identified, for example, by using an algorithm such as the one s« forth in Section 8.2. infrc, 
or another objective algoritiun that objectively ranks genes for use in an ECS stgnatu^. The 
table set forth at Section 8.5 below, which lists signature genes according to titetr ECS 
algorithm score, is therefor particularly useful tor selecting useful and prefe^ed genes m an 
ECS gene signatitre. Hence, as an example, prefetr^l gene sgiantitres of tttis invention 
; comprise nucleic acids that correspond to and/or are capable of detecting at least the firs. 5. 
,0 ,5 20. 25. 30, 35, 40, 45, 50. 60. 70. 80, 90, 100. .120. 120. 130. 140. 150 or all of the 
signature genes Ust in Sections 8.5 - ...... nucleic acids corresponding to and/or capable o 

detecting at least SEQmN0S:l-5; 1-10; M5; 1-20; 1-25; 1-30; 1-35; MO; 1-45; 1- 
50; 1-60; 1-70; 1-80; 1-90; 1-100; l-UO; 1-120; 1-130; 1-140; 1-150; or 1-152. 

^ . . m^/Pntinn mav comprise nucleic acids that 

0 Alternatively, a gene signaiuiv • - 

con«pond «> and/or are capable of detecting orthologs ^d/or homologs otthese genes. A, 
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noted above, homologos and/or orthologs of the rat genes listed in the appendices, infra, can 
be readily identified, for example by their level of sequence identity to the recited rat nucleic 
acid sequences, or by the level of sequence identity and/or homology to amino acid 
sequences that they encode. Alternatively, homologs and orthologs (including those from 
5 other species, such as humans and mice) can be identified by hybridization under conditions 
of appropriate stringency, such as those described, supra, in Section 6.1. In a preferred 
embodiment, appropriate homologs and/or orthologs are identified by routine clustering 
algorithms. For example, the NCBI UniGene database (accessed via the URL 
<http://www.ncbi.nlm.nih.govAJniGene>) groups genes into appropriate clusters of 

10 homologous sequences from the same and/or different species of organism. See, Schuler, 

Mol Med, 1997, 75 (10): 694-698; Schuler et al. Science 1996, 274:540-546; and Boyuski & 
^c\i\x\Qr, Nature Genetics 1995, 10:369-371. 

As a particular example, and not by way of limitation, the gene signature of the 
invention can be used, e.g., in detection kits (for example in any of the applications described 

15 in Sections 7.3-7.5, below) for detecting one or more signature genes. Examples of such kits 
include, but are not limited to, microarrays and PCR reaction kits to name a few. Preferred 
kits will comprise probes and/or primers for amplifying and/or detecting one or more ECS 
signature genes. Hence, a typical kit of the invention will comprises nucleic acids 
(preferably oligonucleotides) that specifically hybridize to ECS signature genes of the 

20 invention under defined conditions (which will depend on the particular application for 
which the kit is used, and include hybridization conditions described in this specification 
supra). 

Kits of the invention will typically comprise at least one probe or primer which is 
capable of specifically hybridizing (e.g., understringent conditions) to an ECS signature 

25 gene. Preferred kits comprises a plurality of probes or primers, each capable of specifically 
hybridizing to a different ECS signature gene. PCR reaction kits will typically comprise a at 
least one pair of primer (z.e., a primer pair) that is capable of specifically amplifying an ECS 
signature gene. Hence, a primer pair in such kits will typically include a "forward" primer 
and a "reverse" primer. Preferred PCR kits of the invention comprise a plurality of primer 

30 pairs, each of which is capable of amplifying a particular ECS signature gene by PCP.. PCR 
kits of the invention also include reverse transcription PCR (RT-PCR) kits, that are capable 
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of converting mRNA from ECS signature genes into cDNA, and then amplifying said cDNA 
in a PGR reaction. Hence, RT-PCR kits of the invention will also typically comprise at least 
one primer pair, and may also comprise a "reverse primer;" Le,, an oligonucleotide that is 
capable of priming the reverse transcription of mRNA from an ECS signature gene. 
5 Preferred RT-PCR kits comprise a plurahty of primer pairs and reverse primers, which are 
capable of reverse transcribing and amplifying mRNAs for a plurality of different ECS 
signature genes. 

hi some embodiments, the reverse primer is a separate and distinct oligonucleotide 
that is different from either the forward or reverse primer. More preferably, however, either 
10 the forward or reverse primer for amplifying an ECS signature gene is also capable of reverse 
transcribing mRNA of that ECS signature gene, and so also ftinctions as a reverse 
transcription primer in an RT-PCR kit. 

Kits may contain additional reagents or other components for detecting ECS signature 
genes, or such additional components may be provided by a user. For example, preferred kits 
15 of the invention may contain at least a polymerase enzyme and one or more labeled 
nucleotides. A kit of the invention may also contain buffer or other reagents. 

7.3. Evaluating Patient Treatment Using ECS Gene Signatures 

The ECS gene signatures identified herein {e.g., in Sections 8.3-8.5, infra) can be 

20 used in evaluating the efficacy of a treatment of a patient suffering from depression, hi this 
evaluation, a peripheral blood sample is taken from the patient before the treatment, and at 1 
day or 1 week after initiating the treatment. Lymphocytes or other accessible peripheral cell 
types are isolated, and the level of expression of all genes within these cells, or a subset of 
genes indicated in Tables 5-7, is then determined by one or more quantitative RNA 

25 measuring techniques. For example, the sample can be treated with Trizol reagent and other 
chemical procedures in order to prepare RNA. Probes and microarrays for the analysis are 
prepared using methods established for gene expression analysis, such as by any one of the 
commonly available quantitative mRNA detection methods, e.g., reverse transcriptase- 
polymerase chain reaction (RT-PCR) or microarray analysis using commercially available 

30 chips. The alteration of gene expression is quantified in the lymphocyte sample(s) and 
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determined to be either higher or lower than the value obtained in the pre-treatment, baseline 
sample(s). 

Resulting data showing gene changes in the same direction as those obtained in the 
ECS samples sets indicates indicate that the treatment method is effective in treating 
5 depression in the patient, whereas a gene change in the opposite direction indicates that the 
treatment is less effective for this purpose. 

7.4. Use Of ECS Gene Signatures To Identify Therapeutic Compounds 

The gene signatures provided in Sections 8.3 - 8.5, infra, can be used in screening for 

10 antidepressive agents. In this screening method, a set of genes is selected as those which 
change significantly following ECS treatments (see Tables 5-7), and the corresponding 
probes are selected by having sequence identity with a significant region of the selected 
genes. Methods and procedures that are common to those skilled in the art of quantifying 
RNA levels in cells, using, e.g., reverse transcriptase-polymerase chain reaction (RT-PCR), 

15 microarray analysis with commercially, high-density cDNA or oligonucleotide-based chips, 
or miniarrays that incorporate relatively small (e.g., 10-300) gene sequences corresponding to 
the genes listed in Sections 8.3-8.5, are used. Microplates wherein each well contains a fixed 
number of cultured neuronal cells from a neuronal, neuroblastoma, astrocytic, or glial cell 
Une, or a primary neuronal or lymphoblast culture, or another suitable cell culture, are 

20 prepared and used as described in the co-pending U.S. patent apiplication Serial No. 
10/175,523, filed June 18, 2002. 

The assay is conducted by contacting neuronal cells with a number of test agents, 
exposing the cells to the test agents for 1-4 days, and the effect on the expression of the select 
set of genes is evaluated. Resuhing data showing gene changes that are in the same direction 

25 as those identified following ECS exposure to rats, particularly those summarized in Tables 
5-7, for a test agent shows that the agent can be effective in treating or preventing depression. 

7.5. Diagnostic Uses of ECS Gene Signatures 

The gene signatures identified in Section 8.3 - 8.5, infra, can also be used for the 
30 prognosis or diagnosis of a patient suspected to be suffering from, or at risk for developing, a 
depression. 
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In this evaluation, samples of blood or other accessible tissues are obtained before, 
and at 1, 2, 4 and additional weeks after, treatment. Circulating lymphocytes (white blood 
cells) or are obtained from the samples and subjected to gene expression analysis using any 
one of the commonly available quantititaive mRNA detection methods, such as reverse 
5 transcriptase-polymerase chain reaction (RT-PCR), microarray analysis using commercially 
available chips, or other methods. The alteration of gene expression is quantified in post- 
treatment sample(s) and determined to be either higher or lower than the value obtained in 
the pre-treatment, baseline sample(s). 

Resulting data showing changes in gene expression that are in the opposite direction 

10 as gene changes associated with ECS as listed in Tables 5-7 indicates that the patient is 

suffering from depression or is at risk for developing a depression, whereas the change in the 
same direction, or absence of a change, indicates a lower risk for depression, hi response to 
an putative antidepressant agent (to determine prognosis of disease following initial baseline 
assessments), changes in gene expression will have prognostic value as well. Li this 

1 5 example of the embodiment of the invention, resulting data that showed changes in gene 
expression that are in the same direction as gene changes associated with ECS as listed in 
Tables 5-7 indicate that the therapeutic agent is improving the patients state of depression, 
whereas the change in the opposite direction, or absence of a change indicates a lower 
effectiveness of the antidepressant potential of the test agent. 

20 
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8. APPENDIX 



8.1. Electroconvulsive Shock (ECS) Therapy and Gene Regulation 
(Literature Review) 



Author, YR 


level/dose 


Post final 
shock tissue 
collection 


Altered Genes and Brain region 


Kim, et al., Mol Cells 2001, 
12:173-177 


Single shock 


3-24 hours 
(hrs) 


Decreased Inositol 1,4,5- 
triphosphate receptor 1 (InsPSRl) 
mRNA in dentate gyrus and CAl 


Jensen et al., / 
Neuroendocrinal 2001, 
13:887-893. 


14 days 




Decreased POMC mRNA in the 
Hypothalamic-pitutitary-adrenal 
(HPA) axis 


Jang et al.. Prog 
Neuropsychopharmacol Biol 
Psychiatry 2001, 25:1571-81 


Chronic 




Adenylate cyclase (AC) isoforms 
in hippocampus and cerebellum 


Kondratyev et al.. Brain Res 
Mol Brain Res 2001, 91 : 1- 
13. 


Low- 
sensitivity 
(minimal) for 
7 days 




ECS treatment prevented elevation 
of Bcl-XS mRNA after prolonged 
seizure activity in hippocampus 
and rhinal cortex 


Koubietal., Brain Res 
2001, 905:161-170 


Acute shock 


72 hrs 


Decreased TH and TpOH protein 
levels in Locus Cereuleus (LC), 
Central tegmental area and Raphe 
Centralis (TpOH only) mRNA 
changes only in LC; increase 
TpOH protein levels in Frontal 
cortex and increase TH protein in 
hippocampus 


Shen et al.. Psychiatry Clin 
Neurosci 2001, 55:75-77 


Acute / 
chronic 


9-24 hrs / 

lhr-2 

weeks 


Decreased Serotonin transporter 
(5-HTT) mRNA in raphe nucleus 


Burnet etal., Eur J 
Pharmacol 2001, 413:213- 
219 


Single shock / 
5 shocks on 
alternate days 




Increased Substance P (tachykinin 
NK(1)) receptor binding site 
density, but no mRNA changes in 
cerebral cortex 


Chen etal.. Synapse 2001, 
39:42-50. 


Acute 
Chronic 


2 hrs 


Increased GRFalpha-1 and 
GFRalpha-2 mRNA in dentate 

gyrus 


Cho et al. , Exp Mol Med 
2000, 32:110-114 


Chronic 




Differential changes in PDE 
isoform mRNA in hippocampus 
and striatum 


Lammers et al., Mol 

Psychiatry 2000, 5:348-388 


10 days 




Increased D3 receptor mRNA 
Shell of the nucleus accumbens 


Madsenetal., Neuroscience 
2000, 98:33-39 


14 days 


24 hrs 


Increased Neuropeptide Y receptor 
Yl and Y2 mRNA in CAl and 
CAB 
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Author, YR 


ECS 
level/dose 


Post final 
shock tissue 
collection 


Altered Genes and Brain region 


Husum et al., 
Neuropharmacology 2000, 
39:1463-73 


lU days 




Increased preproNPY mRNA m 
the CAl and dentate gyrus 


Valentine et al. , Brain Res 
Mol Brain Res 2000, 
75:337-341 




6 hrs 


Increase Fragile X (fmrl) mRNA 
in dentate gyrus granule cell layer 


Burnet et al. , Neurosci Lett 
1999, 277:79-82 


Single, 
repeated. 


Immediately 
and 3 weeks 
after last 
shock 


Decreased 5-HTlA receptor 
mRNA in CA4 Hippocampus and 
increased in dentate gyrus by 
single or repeated ECS. Repeated 
ECS increased cortical 5-HT2A 
receptor mRNA 


Peietal., Neuroscience 
1999, 90:621-627 


Single shock / 
5 shocks over 
lU days 




Single shock = Transient reduction 
of GIRKl mRNA in dentate gyrus 
Repeated shock = increase in 
GIRK2 and GIRKl mRNA in 
dentate gyrus 


Takahashi et al. , J Neurosci 
1999, 19:610-618 






Increased Phosphodiesterase (PDE) 
4 subtypes PDE4A, PDE4B 
mRNA m prefrontal cortex 


Watkins et al., Brain Res 
Mol Brain Res 1998, 61:108- 
113 


Single and 
repeated 


Levels 
return to 
control 
values after 
4onrs 


Dentate gyrus: increased mRNA 
for NMDA subunits NR2A and 
NR2B; reduced mRNA for 
metabotropic glutamate receptor 
mGLUSb 


Suda et al. , J Neurochem 
1998, 71:1554-1563 


Acute and 
repeated 




Increase in Phosphodiesterase type 
IV (PDE IV) mRNA in various 
brain regions after acute shock; 
repeated shock increased PDE IV 
in frontal cortex 


Zetterstrom et al. , Brain Res 
Mol Brain Res 1998, 57:106- 
110 


Acute and 
chronic 


6 -48 hrs 


BDNF increased with both acute 
and chronic, both evident at 6 
hours, chronic shock effete evident 
at 48 hours, but acute shock 
changes gone at 24 hr 


Garcia-Garcia et al. 
Neuroreport 1998, 9:73-77 


Acute and 
chronic 




Corticotropin releasing factor 
(CRF ) increased in paraventricular 
nucleus; propiomelanocortin 
(POMC) increased in arcuate 
nucleus and ventromedial nucleus 
and Proenkephalin (PENK) in 
nucleus accumbens 


Peietal., Neuroreport 
1998, 9:391-394 


5 times over 
10 days 


6 - 24 hrs 


Increased mRNA MAP2 in dentate 
gyrus 
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Author, YR 


ECS 
level/dose 


P^cf finol 

A uSi imsii 
shock tissue 
collection 


Altered Genes and Brain region 


0 ' Donovan et al., / 
lyeurocttCFn lyyo, /u.iZ4i- 
1248 




.5-1 hr and 

A Vim 

4 nrs 


Increased Erg-1 mRNA in dentate 
granule cells returned to normal by 
4 hr Increased Erg-3 peak at 4 hr 


Xing et al. , Brain Res Mol 


Single 


15min-4 hrs 


Immediate increase in 
NurrlmRNA in dentate gyrus 
granule cells, peak at 1 hr and 
return to normal by 4 hr 


Wang et al. , Neurosci Lett 
1996, 226:151-154 






Increased Preprotachykinin-A 
(PPT-A) mRNA in CSF and 
aqueous humor of eye in RABBIT 


Pei et al., Neuroscience 
1997, 78:343-350 


Acute and 
chronic 


6 hrs, 24 hrs 
and 3 weeks 


Acute: reduced voltage-dependent 
potassium channel subunits Kvl.2 
and Kv4.2 at 6 hrs and normal at 
24 hrs in dentate gyrus 
Chronic: Kvl .2 at 6 hrs and 
normal at 24 hrs in dentate gyrus 
Kv4.2 no change at 6 hr elevated 
at 24 nr 


Zachrisson et al . , Eur J 
Pharmacol 1997,319:191- 






Decrease in preprotachykinin-A 
mRNA in caudate -putamen and 
decrease in tachykinin NKl 
receptor mRNA positive neurons 
in Nucleus Accumbens 


Zachrisson et al.. Brain Res 

KAnI Jirniy% 1?^c 1 QQ7 A^X'^XAH 

ivLoi tsrain i\6s i^y/, H-j.j'f /- 
350 


6 shocks in 2 
weeKS 




Decrease in CCK mRNA in 
caudate-putamen 


Jung et al. , Biol Psychiatry 
1996, 40:503-507 


Single and 20 
days 


30 min - 90 
min 


Varying effects in tetradecanoyl 
phorbol acetate-inducible 
sequences (TIS) 1, 7, 8, 1, 21 in 
varying brain regions 


MCvjowan ei ai . , iLur j 
Pharmacol 1996, 306:249-55 


5 shocks over 
5 days 




G protein subunits Gs alpha 
mRNA decreased in CA3 and 
CAl, G9(o) alpha mRNA 

lTir*i*^OG^/^ m /^^nffiiti^ rrx/mo on/I 

iiicrcdocu in aeniaic gyrus ana kjiz 
alpha mRNA reduced in dentate 
gyrus an CA3 


Woldbye et al. , Brain Res 
1996, 720:111-119 


Single shock 


Various 
times up to 
24 hrs 


Prominent induction of c-fos 
mRNA and FOS in dentate gyrus 



-61- 

{M:\323S\100m087us2\000S8277.DOC HHIiilllimnBiOIDIII } 



Author, YR 


ECS 
level/dose 


Post final 
shock tissue 
collection 


Altered Genes and Brain region 


Porter et al . , Brain Res 
1996, 710:97-102 


Single and 5 
shocks over 
l(J days 




Acute: reduced KAl mRNA in 
CA3a increased KA2 MRNA in 
dentate gyrus. 

Repeated shock: reduced KAl 
mKiNA m CAic and UAJc and 
CA3a-b; KA2 mRNA increased in 
dentate gyrus All changes back to 
normal ai 5 weeKs post snocK 


Nay lor et al. , Brain Res Mol 

Drain AcS lyyo, oJ.D*vy'jDj 


Repeated 

snocK 


24 hrs 


Increased mRNA for GluRl 
suDunii or AMr a receptor 


Fitzgerald et al., J 
Neurochem 1996, 66:429- 
432 


Acute 




Increased ICER and CREM 
isoforms in hippocampus, frontal 
cortex and cerebellum 


Chen, et al., Mol 
Pharmacol 1995, 48:880-889 


Acute and 
Chronic 




FosB and FosB-like proteins 


Nibuya et al, J Neurosci 


Acute (Id) 
and Chronic 
(21d) 




Acute: increased BDNF mRNA 
recovered by 18hr Chronic: 
decreased induction of BDNF and 
trkB mRNA, but prolonged their 
expression in dentate gyrus; in 
CAl and CA3, elevated induction 
of BDNF and trkB Mma 


Smith et al., 

Psychopharmacology 1995, 
ion- '111 '5/in 


Acute & 
repeated (5 
shocks in 10 
d) 


4 hrs and 24 
hrs 


Acute and repeated: increased Dl 
and D2 receptor mRNA at 4 hours 
but JNUl 24 hours m Nucleus 
Accumbens 


Burnet et al. , Neuroreport 


5 shocks in 10 
days 




Increased 5-HTlA and 5-HT2A 
receptor mRNA in dentate gyrus 
and neocortex and decreaed 5- 
HTIA in CA3c 


i^inaeiors ei ai . , 
Neuroscience 1995, 65:661- 
670 




1, 3, 9, 24 
hrs 


Brain derived neurotrophic factor 
mRNA and trkB mRNA increased 
in hippocampus 


Dziedzicka-Wasylewska & 

ixu^u^, •/ iydirui 1 rurioffi 

Gen Sect 1995, 102:221-228 






Increased proenkephalin mRNA in 
nucleus accumoens 


Kim et al. , 7 Neurochem 
1994, 63:1991-1994 


Single 


6 hrs & 12 
hrs 


6 hrs: decreased Inositol 1,4,5- 
triphosphate (InsP3) 3-kinase 
mRNA in dentate gyrus returned to 
normal at 12 hours 
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Author, YR 


ECS 
level/dose 


Post final 
shock tissue 
collection 


Altered Genes and Brain region 


Brady et al. , J Qin Invest 
1994, 94:1263-1268 


Single and 
repeated 


24 hrs - 8 
wk 


Single: Increased Corticotrop in- 
releasing hormone (CRH) mRNA 
in the paraventricular nucleus and 
TH mRNA in the locus coeruleus 
at 24 hrs 

Repeated: some increase that 
persisted for up to 3 weeks (CRV) 
and up to o weeks ( 1 ri) 


Passarelli et al., Neurosci 
Lett 1994, 177:147-150 


Single and 
repeated (7d) 


2 hrs 


Increased expression of hcs73 in 
dentate gyrus, CA3 and median 
habenular nucleus; induced hsp70 
mRNA in dentate gyrus (same in 
single and repeated shock) 


Mikkelsen et al., Brain Res 
Mol Brain Res 1994, 23:317- 
322 


14 days 


24 hrs 


lOx increase in PreproNPY mRNA 
in dentate gyrus and piriform 
cortex 


FoUesa et al., Exp Neurol 
1994, 127:37-44 


Corneal shock 
Minimal 
(U.2S, DU- 
70niA 3x in 
Ih) and 
Maximal 
(0.2s, 150 
mA 3x in Ih) 




Min: bFGF mRNA increased in 
entorhinal cortex and hippocampus 
Max: mcreased m bruF mRNA 
in entorhinal cortex, hippocampus, 
olfactory bulb, striatum and 
cerebellum, increased NGF mRNA 
in entorhinal cortex and 
hippocampus 


Kragh et al. , Exp Brain Res 
1994, 98:305-13 


10 shocks, 20 
shocks, or 36 
snocKs given 
on weekdays 
only 


Id, 2d, 30 d 


Increased SS-like and NYP-like 
immunoreactivity 1 and 2 d after 
last shock in outer part of dentate 
molecular layer Most pronounced 
in 36 shock group 


DUUCi Cl al. , J ly6UrOCrl€iTl 

1993, 61:1270-1276 


Acuie ^1 Q or 
3 d) Chronic 
(10dorl4d) 




cnronic . mcreased j-H 1 2 
receptor mRNA 


r adddiciii dnu vjrzi lyciirosci 
Lett 1993, 153:197-201 


Kepeaieu v / u; 
and single 




Repeated: Increased expression of 
somatostantin mRNA in 
hippocampus 


Hosoda and Duman, J 
Neurochem 1993, 60:1335- 
1343 


lOd 




Down-regulation of beta 1 AR 
mRNA in frontal cortex 


K^Du et al Brain Rps: Mnl 
Brain Res 1993, 18:121-126 






illwiCasCU in niI\.l>l/\ aOQ rNr^X 

mRNA in Locus Coeruleus 


Pratt et al. , Brain Res Bull 
1993, 30:691-693 




4- 24 hrs 


Increased GABAa receptor 
subunits alpha 1 and beta 2 mRNA 
in cerebellum at 4-8 hours, 
returned to normal at 24 hrs 
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Author, YR 


ECS 
level/dose 


Post final 
shock tissue 
coucCiion 


Altered Genes and Brain region 


Wong et al. , Brain Res Mol 
Brain Res 1992, 13:19-25 






Induction of the constitutive hsp72 
gene in dentate gyrus 


Lindefors et al., 
Neuroscience 1991, 45:73-80 


5 shocks in 10 
days 




Increased preprocholecystokinin 
and preprotachykinin-A mRNA in 
Edinger-Westphal nucleus 


Kang et al., 

Psychopharmacol Bull 1991, 
27:359-363 


Smgle 


2-8 hrs 


Increased Alpha 1 and gamma2 
GABA receptor subunit mRNA in 
cerebellum and hippocampus not at 
2 hrs, increased at 4 hours 
returned to normal at 8 hours in 
hippocampus 


Cole et al. , . J Neurochem 
1990, 55:1920-1927 


Single and 
repeated 


15 min. and 
4 hrs 


15 min: increased zif/268, c-fos, 
c-jun and jun-B in hippocampus 
and neocortex and pyriform cortex; 
all except c-jun returned to normal 
by 4 hours 


Herman, et al. , Brain Res 
1989, 501:235-246 


7 days 




Increased CRF mRNA in medial 
parvocellular paraventricular 
nucleus (PVN) 


Xeietal., Brain Res Mol 
Brain Res 1989, 6:11-19 


Single and 
repeated (1, 
3, 6 days) 


0.5, 2, 6, 12 
hrs (single); 
24 hrs 
(repeated) 


Single: decrease in DYN mRNA 
at 0.5 hour followed by increase at 
6 hrs in hippocampus; EK mRNA 
increased at 0.5 hour and still high 
at 12 hrs in entorhinal cortex and 
hippocampus 

Repeated: decreased DYN mRNA 
in hippocampus at 24 hrs elevated 
DYN mRNA in striatum and 
hypothalamus. Increased EK 
mRNA in entorhinal cortex at 6 
hrs 


Yoshikawa et al. , ProcNatl 
AdadSci USA 1985, 82:589- 
593 


lOd 




Increased preproenkephalin mRNA 
in hypothalamus and striatum 
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8.2. ECS Gene Selection Algorithm 

Three considerations were quantified and summed to prioritize the significance of 
each gene "hit" as defined in the statistical analysis portion of Section 7.1.1, infra, 
I. Statistical considerations: 

(High score = 27). Score = p value Score + (0.5) Ratio Score + (0.1) Abundance Score 

The significance (p value), fold change (ratio), and gene abundance measured in optical 
intensity units were assigned an integer value as follows: 

p value Score = - log,Q(p value) 

Ratio Score = Ratio if the ratio > 1.5 or 1/Ratio if ratio < 1/1.5 
Abundance Score = log|o (Abundance / 700) 
IL Biological considerations: 

Score = A (10) + B (8) + C (6) + D (5) + E (2) + F (4) (Max. score = 35) 

A. Gene change is in the opposite direction of mRNA or protein change in human 
depression or bipolar disorder, or animal model reported in the literature. 

B. Same as A, but gene change is in the same direction. 

C. Gene is known to be changed by effective treatments for human depression or bipolar 

disorder 

D. Gene is part of a genetic or biochemical pathway known to be associated with human 
depression, bipolar disorder, or their treatment 

E. Gene has been reported to change following seizures in human or animal studies. 

F. Gene's human homologue is in a chromosomal hot spot for depression or bipolar disorder 
as identified by linkage analysis. 

III. Experimental considerations: 

The experimental results score is based on changes in either or both tissues and after acute, 
chronic, or both ECS exposures (Max. score = 14) 

14 Both tissues after acute and chronic 

10 Two tissues, chronic and one tissue acute 

8 Both tissues, chronic 

6 Two tissues acute, one tissue chronic 

4 One tissue chronic, one tissue acute 

2 Two tissues, acute 

1 One tisue chronic 

ECS Algorithm Score ^ 2.2(1) + II + III 

Maximum score possible = 2.2(27)+ 35 + 14 = 108 
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8.3. Fold change (RATIO), significance (p VALUER and ranking of genes 

identified as altered in hippocampus (H) and fi-pntal cortex (F) after either 
actue (A) or chronic fC) ECS exposure 



Gene Name 


Accession 
No. 


RATIO 


P 

VALUE 


Region 

with 
biggest 
change 


Other 
regions 
with 
significant 

change 


Algorithm 
Score 


* cyclooxyg6HciS6 isofoTiti 
(COX-2) 


OO / /ZZ 


1 ,55 




rlA 


A 11 

All 


73.9 


putative potassium chamiel TWIK 


AF022819 


1.83 


9E-04 


HC 


HA 


14.5 


Vascular endothelial growth factor 
(VEGF) 


AA850734 


1.73 


8E-05 


HA 


FA 


11.9 


* brairt'derived neurotrophic factor 
(BDNF) (exon IV) 


D10938 


4,78 


lE-03 


HC 


All 


74.8 


metallothionein 1 A 


AI102562 


2.04 


2E-05 


HA 


FA, HC 


62.1 


Phosphodiesterase 4B, cAMP- 
Specific 


M25350 


1.54 


3E-05 


FA 


FC 


29.8 


fos-related antippti 




3.16 


8E-05 


HA 






neuropeptide Y (NPY) 


M15880 


2.14 


lE-02 


HC 


FC 


54.0 


HMG-CoA reductase 


X55286 


1.48 


3E-05 


HC 


HA, HC 


34.0 


TGFB inducible early growth 
response (TIEGl) 


All /z4 /o 


A /CO 

4.00 


1 t7 AT 

lb-03 


HC 


A 11 

All 


49.0 


nuclear orphan receptor HZF'3 
(NURRl) 


U01146 


1.70 


lE-05 


HA 


HC 


39.9 


JunB 


X54686 


2,21 


lE-07 


FA 


HA 


59.3 


Jun 


AI175959 


1.84 


9E-05 


HA 


FA, FC 


58.8 


* NMDARl glutamate receptor 
subunit 


U11418 


1.59 


2E-02 


HA 




42.0 


* Tissue inhibitor of 
metalloprotemase 1 (TIMPl) 


AI169327 


11.39 


2E-06 


HA 


All 


56,9 


* Narp 


S82649 


5.89 


7E-09 


HC 


All 


46 4 


cytosolic phospholipase A2, group 
IVA 


U38376 


1.79 


6E-07 


HA 


HC 


4!^ 6 


silencer factor B (SF-B) 


X60769 


1.76 


2E-04 


HA 


All 


37.5 


Minoxidil sulfotransferase 


L19998 


0.504 


4E-03 




HA FA 




Vflp nerve prowth factor infiiiriHlp 








HA 

r\j\ 


All 
All 


*fO.D 


phosphatidylinositol 3-kinase p85 
alpha subunit 


D64045 


2.02 


4E-04 


FA 




19.5 


c-fos 


X06769 


2.22 


lE-05 


FC 


HA, HC 


44.9 


Catalase 


AA926149 


0.60 


8E-04 


HA 




44.1 


ceruloplasmin (ferroxidase) 


L33869 


1.81 


2E-05 


FA 


HA, HC 


33.4 


* nerve growth factor-induced 
clone C 


M92433 


2.41 


4E-07 


HA 


HC, FA 


42.5 


thyrotropin-releasmg hormone 


M23643 


2.01 


lE-02 


HC 




42.4 


glial fibrillary acidic protein delta; 
glial fibrillary acidic protein alpha 


AF028784 


1.769 


2E-04 


FA 


HA, HC 


42.0 
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Gene Name 


Accession 
No. 


RATIO 


P 

VALUE 


Region 

with 
biggest 
change 


Other 
regions 
with 
signiflcant 
chance 


Algorithm 
Score 


HES-1 factor 


D13417 


2.09 


5E-03 


HC 


All 


26.7 


activity and neurotransmitter- 
induced early gene protein 4 (ania- 
4) 


AF030089 


2.58 


6E-07 


HA 


FC,HC 


23.8 


Cytochrome P450 Ibl 


AI176856 


2.247 


9E-04 


FA 


HA, HC 


25.3 


arginine vasopressin (Diabetes 
insipidus) 


M25646 


1.52 


4E-03 


HA 


FA 


39.0 


enhancer-of-split and hairy-related 
protein 1 (SHARP 1) 


AF009329 


2.06 


2E-03 


HC 


All 


32.5 


(Ptpnl, non-receptor type 1) 


ATI SOI J.^ 


1 .00 








Zo.j 


Glutathione S-transferase, theta 2 


AI138143 


0.63 


lE-02 


HA 


FA 


26.5 


piuicin piiubpiiaiaoc 1, icguiaiory 
(inhibitor) subunit 1 A 




1 A'> 

i.OZ 






rA, ilA 


22.7 


krox20; (early growth response 2) 


U78102 


3.13 


7E-04 


FA 


HA, HC 


36.2 


iti^'Aflf^iiVriTi ft YP^r^fi^niTiv 
ililvXlCUJviil U ICWCpiUl 


iViJO JO / 






P A 

rA 




22.2 


S-adenosylmethionine 
decarboxylase 


AI008131 


1.55 


lE-02 


HA 


FA 


24.8 


ets variant gene 3 


AA859750 


1.62 


lE-05 


HA 


All 


28.5 


* neiu^l receptor ^BiiPy^gllF™ 
kinase (TrkB) 


M55291 


1.52 


2E-02 


HA 


HC, FA 


20.4 


*Vesl 


AB003726 


6.50 


2E-07 


HA 


FA, HC 




pyruvate dehydrogenase 
phosphatase isoenzyme 1 


AF062740 


1.83 


2E-02 


HC 


HA, FA 


17.6 


F*^1T1 QSQS7 {M^r\A Vinmnlna 

(MAD4)) 




u.o / 


1 Th CiA 


rlA 


A 11 

All 


27.2 






1 .o\)o 




tf A 




M.I 


protein 2 (GRB2) 


AT 17^77 A 


1 .06 


Ap 


T-I A 

rlA 


ri\^ 


32. U 


prostaglandin D synthase 


J04488 


0.35 


6E-03 


HC 


PC 


25.5 


inciduuLropic giuiamaic rccepior 
(GluRl) 


ivioiuyy 


u.ou 




XJ A 

rlA 




25.5 


Decorin 


AI639233 


0.66 


2E-02 


HA 




31.5 


aldehyde dehydrogenase 


M73714 


0.45 


2E-05 


HA 




20.6 


rxuuyuys ivanus norvegicus ciyiNA 
clone 




U.OZ 


ITS 


HA 


HC, FA 


17.3 


activin type I receptor 


L19341 


1.63 


3E-05 


HC 




7.2 


growth factor (Arc) 


U19866 


4.82 


6E-09 


HA 


FA,HC 


29.0 


Interferon-related developmental 
regulator 1 


AI014163 


2.03 


6E-03 


HC 


HA, FA 


21.7 


neuritin 


U88958 


1.56 


lE-08 


HA 


FA,HC 


27.7 


myr 6 myosin heavy chain 


U60416 


0.65 


3E-05 


HA 


HC 


26.7 


interferon gamma receptor (Ifiigr) 


U68272 


0.65 


5E-04 


PC 




20.1 
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Gene Name 


Accession 
N . 


RATIO 


P 

VALUE 


Region 

with 
biggest 
change 


Other 
regions 
with 
significant 
change 


Algorithm 
Score 


lost on transformation 1 (LOTl) 


U72620 


1.58 


8E-05 


FA 




6.0 


DNAJ (Hsp40) homolog, subfamily 
B, member 5 


AA891542 


2.32 


1E.07 


HA 


HC, FA 


25.7 


JE product (small inducibe 
cytokine A2) 


X17053 


1.85 


3E-03 


FA 


HC, HA 


25.4 


Activity and Neurotransmitter 
induced early gene 3; (ania-3) 


AF030088 


4,787 


4E-05 


HC 


HA 


24.9 


Roaz 


U92564 


0.65 


5E-05 


HA 


HC 


14.3 


3LHi 54/LLluu FlFase (protein 
tyrosine phosphatase, non-receptor 
16) 


CO 1 A^O 

S81478 


1.67 


8E-03 


tT A 

HA 


FA 


24.3 


matrilysin ( Matrix 
metalloproteinase 7 ) 


L24374 


1.62 


9E-03 


HA 




4.3 


EST194857 


AA891054 


1.82 


5E-03 


FC 


All 


24.0 


tachykmm 2 (Neurokinin B 
precursor (NKB)) 


A>ri 1 n 
M1041U 










22.1 


rx04104s (Neurodl) 


AI639109 


0.42 


2E-03 


HA 


FA, HC 


21,8 


EST190175 (atrophin-1 related 
protein) 


AA800678 


0.64 


9E-06 


HA 


HC, FA 


21.5 


sniiitp rarripr familv mpmbpr 1 




0 48 






FP HA 


1 1 A 


■flaviti-cnntflifiiniT mntinnwiypna^ip 1 

J-IO-V 11.1 wlJAiUCI nil life llldlll^V/A Y^WliaOv' 1. 


M84719 


0 67 




FA 


FT A 


1 0 1 


SPA-1 like protein 


Mlil516 


1.57 


lE-02 


HA 




42 


PPlMMllO 


S74907 


1.76 


3E-03 


HC 




19.9 


leucine-rich acidic nuclear protein 




1.60 


at: 

4Ji-02 


FC 




5.0 


BHF-1 


D82074 


0.44 


6E-04 


HA 


FA. HC 


19.0 


rx01019s Rattus norvegicus cDNA 
clone 


AI639256 


0.62 


4E-04 


HA 


HC 


12.8 


Nuclear receptor subfamily 4, 

group A, member 3 


AI176710 


1.62 


8E-04 


FA 


HA 


10.8 


Isovaleryl Coenzyme A 
dehydrogenase 


All 02838 


0.64 


3E-03 


HC 




18.4 


insulin-like growth factor II 
(somatomedin A) 


X17012 


0.44 


5E-02 


HC 


FC 


17.8 




UiOloo 


3.16 


OT7 AC 

8E-05 


HA 




17.0 


UI-R-EO-bs-h-03-O-UI.sl UI-R-EO 
Rattus norvegicus cDNA clone 


AA859627 


1.56 


2E-02 


HA 




5.0 


Calvasculin; SI 00 calcium-binding 
protein A4 


X06916 


1.58 


3E-02 


HC 




16.8 


EST189184 


AA799687 


0.65 


2E-03 


HA 




4.8 


developmentally regulated protein 
TPOl 


AIO 12275 


0.61 


4E-04 


HA 


HC 


10.6 


UI.R-AO-bd-e-03-O-UI.sl UI-R-AO 
Rattus norvegicus cDNA clone 


AA866485 


1.53 


lE-03 


FA 


HA 


10.0 
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Gene Name 


Accessi n 
No. 


RATIO 


p 

* 

VALUE 


Region 

with 
biggest 


Other 
regions 
with 

cicvnifiriiTif 

change 


Algorithm 

Score 


cyclic AMP-regulated 
phosphoprotein 


bOjUV 1 


z.Uy 


'Trr A*? 

/b-U/ 


TLT A 
HA 


HC 


15.8 


di-N-acetylchitobiase 


My J /oo 


U.OO 


ib-U5 


TJ" A 

HA 




5.8 


mothers against dpp 3 homolog 


Uoo4 /y 


1 TO 

1. /o 


3b-u3 


TJ A 

HA 




15.7 


lib Izz^z4i (blUU calcium bmaing 
protein A 1) 


Ai2zoj4o 


0.65 


T "C? AO 

3b-02 


HC 




15.7 


UCAAi/eniiancerouiaing, protein 
(C/EBP) delta 


Mo5149 


2.95 


2E-05 


FA 


FC 


15.4 


lio 1 1 yoo 14 (tescaicinj 




1 CI 

1.53 


1 "C A'> 


HC 


tl A 

HA 


9.1 


EST220045 


All 76460 


7.81 


3E-01 


FC 




15.0 


Potassium voltage gated channel, 
Shai-related family member 2 


AI230211 


2.437 


5E-02 


HA 


FC 


14.3 


transducin-like enhancer of split 
l(Tlel) 


A A QKf\QA 


U.oU 


1 A/C 


HA 


"E7 A 

FA 


14.0 


♦Fibroblast growth factor receptor 
1 


D 12498 


1.51 


2E-03 


HA 




7.7 


chemokine CX3C 


AFO30358 


1.59 


lE-05 


HA 


HC 


12.9 


neuropilin 


AF016296 


1.67 


2E-03 


HA 


HC 


12.8 




AAoyizUo 


1 iO 

1.63 


A AC 

4E-05 


HA 


HC 


12.8 


Xio 1 IVUjjU 


A A CHAQ^I 


1.55 


T i:; AC 


r? A 


HA 


12.4 


Ho I z 1^103 ^^rJvDUO oinamg 
protein 3) 


ATI f\1QnA 

AllUio /4 


U. /z 


1"C A/i 


HC 


TT A 

HA 


12.0 


DNA-damage-inducible transcript 
1 


L32591 


1.68 


2E-03 


HC 




5.8 


thymosin beta-10 


M58404 


1.51 


3E-04 


HC 




3.8 ^ 


Y box protein 1 


AI230572 


1.53 


7E-03 


HA 




3.6 


phosducin-like protein 


L15354 


0.56 


6E-06 


HA 


FA 


11.3 


Ash-s 


D49847 


1.56 


lE-03 


HC 


HA 


11.2 


connexin protein Cx26 (gap 
junction beta-2 protein) 


1 1 < 


U.OO 


1 t? A*! 


T? A 

FA 


HC 


11.2 


Hariy growui response i ^liijKi ^ 


ATI 'TiC/C/CO 
All /OOOZ 


l.oi 


1 "C A1 

lb-03 


U A 

HA 


FA 


11.1 


Jboizijujj (suxuiarto vj33; 


A T1 /ic\nc/z 

Alloy/jo 


1.78 


3E-03 


FA 


HA 


10.7 


rx04422s Rattus norvegicus cDNA 
clone 


AI639169 


0.66 


5E-04 


HA 




4.6 


protein tyrosine phosphatase 
(PtP4a2) 


AJ007016 


1.51 


3E-03 


HA 


FA 


10.4 




A AXOl 7^^^ 
/VrVOy 1 / J J 


1. JO 


Jli-U4 


rlA 






rx04826s Rattus norvegicus cDNA 
clone (similar to Nedd4) 


AI639058 


1.62 


2E-02 


HC 


HA 


10.0 


UI-R-E0-€i-e-02-0-in.sl UI-R-EO 
Rattus norvegicus cDNA clone 


AA874943 


0.64 


3E-04 


HA 


HC 


9.2 


Transcription factor UBF 


M61725 


1.58 


8E-03 


HA 


HC 


8.8 
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Gene Name 


Accession 
No. 


RATIO 


p 

VALUE 


Regi n 

with 
biggest 
change 


Other 
regions 
with 
significant 
change 


Algorithm 
Score 


EST203S92 (cnromobox bomolog 
1) 


AI009141 


0.63 


3E-02 


HC 


HA 


8.7 


Lysyl oxidase 


S77494 


2.407 


6E-03 


HC 




8.5 


Small proline-rich protein gene 
(comifin alpha) 


AA891911 


1.52 


2E-02 


HC 




6.9 


EST196165 


AA892362 


0.65 


3E-04 


HA 




6 7 


ESTl 06597 


H31990 


0.64 


lE-03 


HC 




6.6 


heme oxygenase-3 


AF058787 


1.61 


9E-03 


HC 




6.5 


CO 1 ly \jsJ\\j^L nomoiogj 




U.04 


yb-U4 


riA 




6.3 


protein kinase MNK2 


Z21935 


1.81 


2E-02 


HA 




6.2 


SHB (Src homology 2 domain 
containing) adaptor protein B 


AA859468 


1.77 


7E-03 


HC 


HA 


5.2 


rx00382s Rattus norvegicus cDNA 
clone 


A TAIOl 
AlOJV 1 JJ 


U.JO 




tJ A 

rlA 




5.2 


UI-R-EO-bx^e-ll-0<Ul.sl UI-R-EO 
Rattus norvegicus cDNA clone 


AA859690 


1.94 


3E-08 


HC 




4.7 


EST190198 (coagulation factor III) 


AA800701 


1.58 


2E-04 


HA 




4.4 


EST190228 


AA800731 


0.61 


7E-03 


FC 




4.1 


^^181^197395 (reticulocalbin,2). ^ 


AA893592 


1.67 


5E-02 


HC 




4.1 


nuclear pore complex protein 
NupSO 


U41845 


1.52 


2E-04 


HA 




3.9 


♦clathrin, light polypeptide (LCB) 


AA874955 


1.62 


4E-03 


HA 




3.8 


XT A 1 t^Vtf^-t-a-i-n 

fNAv-^-i proxein 


A "cni <Q1 1 


1 HI 

1.03 




U A 

rIA 




3.8 


Dcia ucicnsm*! 


ArUOooOU 


1 Oil 


ZH-UZ 


TJ A 

rlA 




5.1 


voltage -dependent potassium 
chaimel; (shaker-related subfamily) 




1 CI 

l.oi 




TT A 

rlA 




3.5 


rx00967s Rattus norvecficus cDNA 
clone 


AI639015 


0.61 


4E-04 


HA 




'\ 1 

J. 1 


Beta-carotene 15, 15'-dioxygenase 


AI014135 


0.59 


2E-02 


FA 




2.9 


BRL-3 A binding protein 


A09811 


0.32 


lE-02 


HC 




2.9 


potassitmi voltage gated channel, 
shaker related subfamily, member 4 


M32867 


1.61 


lE-02 


HC 




2.7 



-70- 

{M:\3235\l 00m087us2\000S8277.DOC IHHIIDiEiniiMIMDIIillllD } 



8.4. Supplemental list of highly specific genes that change less than 1.5-fold 
(Columns are as defined in Section 8.3. supra) 



Gene 


Accession 
No. 


RATIO 


P 

VALUE 


Region 
with 

lowest p 
value 


Other 
regions 
with 
signiOcant 
change 


Alg rithm 
Score 


Bcl-2 associated death agonist 
(BAD) 


AA8 18072 


0.8 


3E-05 


HA 


HC 


32 


beta-tubulin T beta 15 (aa 1- 
445) 


X03369 


1.27 


lE-05 


HC 




8 


diacylglycerol kinase 


S49760 


137 


4E-04 


HC 




7 


EST 189694 Nomialized rat 
heart 


AA800197 


0.69 


4E-07 


HA 


HC 


11 


EST 195383 Nonnalized rat 

Icidnev 


AA891580 


0.67 


2E-05 


HC 


HA 


9 


EST 197783 Nonnalized rat 
nlacenta 


AA893980 


0.76 


6E-06 


HA 


HC 


10 


GAB A-A receptor 


X51992 


0.79 


2E-04 


HC 


HA 


19 


Glypican 


L02896 


1.28 


3E-03 


HC 


FA, FC 


13 


Growth factor receptor bound 

nrntpin 2 


AA801130 


1.23 


5E-04 


HC 


HA, FC 


16 


Light molecular weight 
neurofilament 


AF031880 


1.4 


lE-05 


HA 


HC 


24 


XT X"l A. A. ' ' 

Neurofilament protein middle 
(NF-M) 


Z12152 


1.26 


4E-08 


HA 




23 


Neurofilament protein; rat 
smallest (NF-L) 


M25638 


1.49 


2E-07 


HA 


HC 


22 


Neuron glucose transporter 


D 13962 


1.4 


7E-05 


FA 




5 


Phospholipase C-beta 1 


L14323 


0.79 


3E-05 


HA 




24 


Phospholipase D gene 2 


NM_033299 


0.72 


lE-04 


HA 


HC 


23 


Ras homolog enriched in brain 


Al 009605 


1.42 


3E-03 


HA 


HC 


12 


SH3-domain GRB2-like 1 


AlOlOllO 


1.26 


1E.05 


HA 




12 


Sodium-dependent 
neurotransmitter transporter, 
ventral midbrain 


S56141 


0.78 


lE-05 


HA 




5 


Synapsin 2 


M27925 


1.34 


4E-04 


FC 




6 
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8.5. Confirmation of Select Gene Changes by Quantitative RT-PCR 
All samples for ECS-treated groups and sham control groups used for microarray 
studies were analyzed for each gene (n = 8-10 per group). 



Frontal Cortex, Acute ECS 







QRT-PCR 


Microarray 


UniGene 


Name 


Ratio 


p value 


Ratio 


p value 


Rn.32777 


Ceruloplasmin (ferroxidase) 


1.35 


2.8E-02 


1.81 


2.0E-05 


Rn.10599 


Phosphoinositide 3-kinase, 
regulatory subunit, 
polypeptide 1 


1.59 


5.8E-04 


2.02 


4.0E-03 


Rn.6977 


Lost on transformation 1 


1.90 


2.1E-04 


1.58 


8.0E-05 


Rn.1716 


Interleukin 6 receptor 


2.24 


3.1E-08 


1.74 


3.0E-05 















Frontal Cortex, Chronic ECS 







QRT-PCR 


Microarray 


UniGene 


Name 


Ratio 


p value 


Ratio 


p value 


Rn.11400 


Prostaglandin D synthase 


0.61 


1.0E-02 


0.71 


0.03 


Rn.7730 


Prostaglandin E synthase 


2.18 


1.5E-02 


not on U34A chip 


Rn.506 


Synapsin 2 


1.29 


1.0E-03 


1.34 


0.00 















Hippocam 


pus, Acute ECS 






QRT-PCR 


Microarray 


UniGene 


Name 


Ratio 


p value 


Ratio 


p value 


Rn.87787 


glutamate receptor, 
metabotropic 1 


0.63 


2.5E-06 


0.60 


3.0E-05 


Rn.19727 


Hairy and enhancer of split 1 
(Drosophila) 


0.71 


3.3E-02 


1.55 


4.0E-05 


Rn.3812 


Normalized rat heart, Bento 
Soares Rattus sp. cDNA 
RHEAE07 


0.72 


1.4E-05 


0.65 


2.0E-03 


Rn.9113 


Aldehyde dehydrogenase 

family 3, subfamily A2 


0.86 


5.6E-03 


0.44 


2.0E-05 


Rn.3360 


Growth factor receptor 
bound protein 2 


1.36 


3.1E-02 


1.68 


6.0E-06 


Rn.10815 


Neuropilin 


1.36 


1 .2E-04 


1.67 


2.0E-03 


Rn.9797 


Fibroblast growth factor 
receptor 1 


1.44 


1 .6E-03 


1.51 


2.0E-03 


Rn.2045 


Protein tyrosine 
phosphatase type IVA, 
member 2 


1.46 


2.0E-03 


1.51 


3.0E-03 


Rn.3546 


Neuritin 


1.59 


4.9E-06 


1.56 


1.0E-08 


Rn.9839 


Nuclear receptor subfamily 
4, group A, member 2 


1.76 


6.8E-05 


1.70 


1.0E-05 


Rn.44320 


Avian sarcoma virus 1 7 (v- 
jun) oncogene homolog 


1.80 


1.7E-07 


1.84 


9.0E-05 


Rn.859 


Ras homolog enriched in 

brain 


1.83 


1.0E-07 


1.42 


1.0E-05 


Rn.11183 


DNA-damage inducible 
transcript 3 


1.89 


1.6E-08 


1.78 


4.0E-03 


Rn.3723 


Interferon-related 


1.93 


1.3E-07 


2.01 


3.0E-05 
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developmental regulator 1 










D82074 


Rattus sp. mRNA for BHF-1 


1.97 


7.0E-04 


0.44 


6.0E-04 


Z21935 


R.norvegicus protein kinase 
rMNK2 


2.28 


1.6E-02 


1.81 


2.0E-02 


S65091 


cyclic AMP-regulated 
phosphoprotein 


2.52 


6.7E-05 


2 09 


7 OE-07 


Rn.2398 


TGFB inducible earlv arowth 
response 


2.79 


1 .4E-08 


3 62 


1 OF-06 

1 .\/L_ \J\J 


Rn 40517 


Activitv and 

neurotransmltter-induced 
early gene protein 4 (ania-4) 


3 64 


2.9E-10 


2 58 


6 OE-07 


Rn.44369 


Prostaglandin-endoperoxide 
synthase 2 


4.03 


1.5E-07 


7.33 


9.0E-08 


Rn.25754 


Tissue inhibitor of 
metalloproteinase 1 


4.71 


8.7E-08 


11.39 


2.0E-06 


Rn.11266 


Brain derived neurotrophic 
factor 


5.73 


2.9E-08 


3.85 


2.0E-04 


S82649 


neuronal activity-regulated 

pentraxin 


7.24 


2.0E-10 


5.89 


7.0E-09 


Rn.37500 


Honner, neuronal immediate 
early gene, 1 


7.42 


2.5E-07 


6.50 


2.0E-07 



Hippocampus, Chronic ECS 







QRT-PCR 


Microarray 


UniGene 


Name 


Ratio 


p value 


Ratio 


p value 


Rn.964 


Insulin-like growth factor II 
(somatomedin A) 


0.33 


4.4E-02 


0.44 


5.0E-02 


Rn.11400 


Prostaglandin D synthase 


0.34 


3.9E-02 


0.35 


6.0E-03 


Rn.19727 


Hairy and enhancer of split 1 
(Drosophila) 


1.27 


1.8E-02 


2.09 


5.0E-03 


M25638 


Rat smallest neurofilament 
protein (NF-L) mRNA 


1.31 


2.1E-03 






Rn.7044 


Glypican 1 


1.33 


1.6E-02 






Rn.78188 


heme oxygenase-3 (HO-3) 
mRNA 


1.42 


2.3E-02 


1.61 


9.0E-03 


Rn.859 


Ras homolog enriched in 

brain 


1.54 


6.7E-04 






Rn.9756 


Protein phosphatase 1 , 
regulatory (inhibitor) subunit 
1A 


1.58 


4.8E-06 


1.62 


6.0E-03 


Rn.3546 


Neuritin 


1.58 


8.6E-05 


1.44 


8.0E-04 


Rn.31799 


Pyruvate dehydrogenase 
phosphatase isoenzyme 1 


1.72 


1.3E-03 


1.83 


2.0E-02 


Rn.9714 


Neuropeptide Y 


2.40 


5.2E-06 


2.14 


1.0E-02 


Rn.22 


Thyrotropin releasing 
hormone 


21.42 


1.8E-03 


2.01 


1.0E-02 
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§ 



g 



u 

'S 

CO 

> 



§ 

o 

s 

-4-* 

CJ 
W 



00 



Algorithm Scores 


ECS 
Algorithm 
Score 


87.9 


76.5 


75.9 


74.8 


68.1 


67.8 


65.9 


CD 


CD 




59.9 


59.7 


59.3 


Groups 1, II 
(Literature) 
and III 


73.9 


14.5 


11.9 


74.8 


62.1 


29.8 


59.9 


CO 


LO 




39.9 


1^ 


59.3 


PGI II 




CVJ 
CO 


CO 


o 


CD 


00 
CO 


CD 


o 

CO 


O 


CVJ 


o 

CSJ 


CVJ 
ID 


o 


PGI Data, Group II: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


o 


o 


CM 


o 


CD 


o 


CD 


o 


O 


o 


CVJ 




o 


Human 
cells 


CD 


CVJ 


CD 


o 


O 


o 


O 


CD 


O 


CVJ 


o 


CVJ 


o 


Human 
tissue 


00 


o 

ID 


CD 


o 


O 


00 
CO 


O 


CSJ 


O 


o 


00 


o 


o 


Unigene number 


Mouse 
Unigene 


CO 
CO 

E 
^ 


Mm.10800 j 


Mm.31540 


Mm. 1442 


Mm.192991 


T— 

00 

o 

CVJ 

E 
^ 


NA 


NA 


Mm. 154796 


Mm.4292 


o 
in 

CO 

E 


h- 
in 

CO 

E 
^ 


NA 


Human 
Unigene 


Hs.1 96384 


Hs.79351 


Hs.73793 


Hs.56023 


Mm.1929 
91 


Hs.188 


NA 


Hs.11899 


Hs.1832 


Hs,82173 


Hs.82120 


Hs.748 


< 


Rat 
unigene 


Rn.44369 


Rn. 15693 


Rn.1923 


Rn. 11266 j 


Rn.54397 


Rn.2485 


Rn.3212 


Rn.9437 


Rn.9714 


Rn.2398 


Rn.88129 


Rn.9797 


< 
z 


Accession 
No, 


S67722 


AF022819 


AA850734 


AI030286 


Ah 02562 


M25350 


AA875032 


X55286 


Ml 5880 


AI071299 


U01146 


D12498 


X54686 


Gene Name 
(SEQ ID NO.) 


cyclooxygenase isofomi 

(COX-2) 

(SEQIDNC:1) 


putative potassium 
channel TWIK 
(SEQ ID N0:2) 


Vascular endothelial 
growth factor (VEGF) 
(SEQ ID N0:3) 


brain-derived neurotrophic 
factor (BDNF) (exon IV) 
(SEQ ID N0:4) 


metallothionein 1 A 
(SEQ ID N0:5) 


Phosphodiesterase 4B, 
cAMP-SpecifIc 
(SEQ ID N0:6) 


fos-related antigen (fra-2) 
(SEQ ID N0:7) 


HMG-CoA reductase 
(SEQ ID N0:8) 


neuropeptide Y (NPY) 
(SEQ ID N0:9) 


TGFB inducible early 
growth response (TIEG1) 
(SEQ IDNO:10) 


nuclear orphan receptor 
HZF-3 (NURR1) 
(SEQ ID N0:11) 


FGF receptor-l 
(SEQIDN0:12) 


Jun B 

(SEQ ID N0:13) 



o 
q 

(N 

00 
irj 
O 

o 

CN 

CO 

r-- 

C30 

1 



CO 
CM 



Algorithm Scores 


ECS 
Algorithm 
Score 


58.8 


00 

in 


56.9 


56.4 


54.6 


49.5 


47.6 


46.5 


45.5 


44.9 


44.1 




43.4 


42.5 


Groups 1, II 
(Literature) 
and III 


58.8 


CM 


56.9 


46.4 


48.6 


37.5 


25.6 


46.5 


19.5 


44.9 


44.1 


CO 


33.4 


42.5 


PGI 11 


o 


CD 


o 


o 


CD 


CM 


CM 
CM 


o 


CD 
CM 


o 


o 


o 


o 


o 


PGI Data, Group II: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


o 


o 


o 


o 


o 


CD 


CM 


o 


O 


o 


o 


o 


o 


o 


Human 
cells 


o 


o 


o 


o 


CD 


CD 


O 


o 


O 


o 


o 


o 


o 


o 


Human 
tissue 


o 


CD 


o 


o 


O 


O 


o 


o 


CD 
CM 


o 


o 


o 


o 


o 


Unigene number 


Mouse 
Unigene 


NA 


Mm.3292 


VN 


NA 


CD 
00 

E 
^ 


Mm.4863 


Mm.17339 


NA 


Mm.3058 


VN 


Mm.4215 


Mm .42083 


Mm.13787 


NA 


Human 
Unigene 


NA 


Hs.105 


NA 


NA 


Hs.211587 


Hs.99029 


Hs.142 


Hs.171014 


Hs.6241 


NA 


Hs.395771 


Hs.41143 


Hs.296634 


NA 


Rat 
unigene 


Rn.44320 


Rn.9840 


Rn.25754 


NA 


Rn.10162 


Rn.6479 


Rn.1507 


Rn.9704 


Rn.10599 


NA 


VN 


Rn.45523 


Rn.32777 


NA 


Accession 
No. 


AM 75959 


U11418 


AM 69327 


882649 


U38376 


X60769 


LI 9998 


M74223 


D64045 


X06769 


AA891848 


L14323 


L33869 


M92433 


Gene Name 
(SEQ ID NO.) 


Jun 

(SEQIDN0:14) 


NMDAR1 glutamate 
receptor subunit 
(SEQIDN0:15) 


Tissue inhibitor of 
metalloproteinase 1 
{TIMP1) 

(SEQIDN0:16) 


Narp 

(SEQIDN0:17) 


cytosolic phospholipase 
A2, group IVA 
(SEQ IDN0:18) 


silencer factor B (SF-B) 
(SEQ ID N0:19) 


Minoxidil sulfotransferase 
(SEQ ID NO:20) 


VGF nerve growth factor 

inducible 

(SEQIDN0:21) 


phosphatidylinositol 3- 
kinase p85 alpha subunit 
(SEQ ID NO:22) 


c-fos 

(SEQ ID NO:23) 


Catalase 
(SEQIDN0:24) 


phospholipase C-beta 1 
(SEQ ID NO:25) 


ceruloplasmin 
(ferroxidase) 
(SEQ ID NO:26) 


nerve grovrth factor- 
Induced clone C 
(SEQ ID NO:27) 
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00 
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CM 

00 

1 
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ro 
CM 
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Algorithm Scores 


ECS 
Algorithm 
Score 


CO 
CO 


in 

CO 


34.8 


34.5 


34.4 


CO 


33.9 


33.6 


33.2 


32.2 


CNJ 
CO 


CN 
CO 


CN 
CO 


Groups 1, II 
(Literature) \ 
and III 1 


o 

CO 


<j> 

CNJ 


24.8 


28.5 


20.4 


00 
CNJ 


33.9 


17.6 


27.2 


32.2 


CNI 
CO 


CN 
CO 


CN 
CO 


PGI II 


CO 


CD 


o 


CD 




CD 


o 


CD 


<D 


o 


o 


O 


O 


PGI Data, Group 11: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


o 


o 


o 


CD 


o 


O 


o 


O 


CD 


o 


o 


O 


O 


Human 
cells 


CO 


CD 


o 


O 


CD 


CD 


o 


O 


O 


o 


o 


O 


O 


Human 
tissue 


o 


o 


o 


O 


00 


O 




CD 


O 


o 


o 


O 


O 


Unigene number 


Mouse 
Unigene 


Mm.1956 


Mm. 142140 


Mm.7880 


NA 


Mm.3993 


Mm.1956 


NA 


NA 


NA 


VN 


Mm.4387 


Mm.6900 


Mm.6900 


Human 
Unigene 


Hs.211584 


Hs.71346 


Hs.262476 


NA 


Hs.47860 


Hs.211584 


NA 


Hs.22265 


NA 


NA 


Hs.76366 


Hs.381152 


Hs.381152 


Rat 
unigene 


Rn, 18568 


Rn.10971 


Rn.29949 


Rn.7937 


Rn.11246 


Rn.18568 


Rn.37500 


Rn.31799 


Rn.3279 


Rn.40130 


Rn.36696 


Rn.3360 


Rn.3360 


Accession 
No. 


AF031880 


Z12152 


AI008131 


AA859750 


M55291 


M25638 


AB003726 


AF062740 


AA892154 


AI235890 


AA818072 


All 70776 


AA801130 


Gene Name 
(SEQ [D NO.) 


light molecular weight 
neurofilament 
(SEQ ID NO:39) 


neurofilament protein 
middle (NF-M) 
(SEQ ID NO:40) 


S-adenosylmethionine 
decarboxylase 
(SEQ ID N0:41) 


ets variant gene 3 
(SEQIDNO:42) 


neural receptor protein- 
tyrosine kinase (TrkB) 
(SEQ ID NO:43) 


neurofilament protein; rat 
smallest (NF--L) 
(SEQ ID NO:44) 


VesI 

(SEQ ID NO:45) 


pyruvate dehydrogenase 
phosphatase isoenzyme 1 
(SEQ ID NO:46) 


EST1 95957 (Mad4 
homolog) 
(SEQ ID NO:47) 


MHC class IbRTI .S3 
(SEQ ID NO:48) 


Bad: bcl-2 associated 
death agonist 
(SEQ ID NO:49) 


Growth factor receptor 
bound protein 2 (GRB2) 
(SEQ ID NO:50) 


Growth factor receptor 
bound protein 2 (GRB2) 
(SEQ ID N0:51) 
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fS 
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3 
00 
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ro 

CN 



Algorithm Scores 


ECS 
Algorithm 
Score 


31.5 


31.5 


31.5 


CO 


30.6 


29.3 


29.2 


a> 

CM 


27.7 


27.7 


26.7 


26.1 


CO 
CM 


CD 
CM 


Groups 1, II 
(Literature) 
and III 


31.5 


25.5 


25.5 


ir> 

CM 


20.6 


17.3 


CM 


<J> 

CM 


21.7 


27.7 


26.7 


20.1 


CD 


CO 


PGI II 


o 


CD 


CD 


CD 


o 


CM 


CM 
CM 


O 


CD 


o 


o 


CD 


O 
CN 


O 

T- 


PGI Data, Group il: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


o 


O 


CD 


o 


o 


CM 


CD 


O 


CO 


o 


o 


O 


O 


o 


Human 
ceils 


o 


CO 


O 


CD 


o 


O 


O 


O 


o 


o 


o 


CD 


CM 


o 


Human 
tissue 


o 


O 


O 


o 


o 


O 


CD 


O 


o 


o 


o 


O 


CO 


o 


Unigene number 


Mouse 
Unigene 


Mm.56769 


Mm.157764 


Mm.1008 


Mm.261561 


Mm.4210 


NA 


Mm.689 


Mm. 25405 


00 
CD 

E 


VN 


Mm.3536 


Mm.549 


Mm.220978 


Mm.20892 


Human 
Unigene 


Hs.433989 


Hs.32945 


Hs.430637 


Hs.24969 


Hs.1 59608 


NA 


Hs.1 50402 


NA 


Hs.7879 


NA 


NA 


Hs.180866 


Hs.75825 


Hs.6439 


Rat 
unigene 


Rn.3819 


Rn.87787 


Rn.11400 


Rn.10368 


CO 

d 


CD 

o> 

ID 
CD 

T- 

d 


Rn.87899 


CD 
CO 
O 
O 

d 


Rn.3723 


Rn.3546 


Rn.10640 


Rn.19927 


Rn.6977 


Rn.506 


Accession 
No. 


AI639233 ^ 


M61099 


J04488 


X51992 


M73714 


AI638960 


LI 9341 


U 19866 


AI014163 


U88958 


U60416 


U68272 


U72620 


M27925 


Gene Name 
(SEQ ID NO.) 


Decorin 

(SEQ ID NO:52) 


metabotropic glutamate 
receptor (GluRI) 
(SEQ ID NO:53) 


prostaglandin D synthase 
(SEQ ID NO:54) 


gamma-aminobutyric acid 
A receptor, alpha 5 
(SEQ ID NO:55) 


aldehyde dehydrogenase 
(SEQ ID NO:56) 


rx00909s Rattus 
norvegicus cDNA clone 
(SEQ ID NO:57) 


activin type 1 receptor 
(SEQ ID NO:58) 


growth factor (Arc) 
(SEQ ID NO:59) 


Interferon-related 
developmental regulator 1 
(SEQ ID NO:60) 


neuritin 

(SEQIDN0:61) 


myr 6 myosin heavy chain 
(SEQ ID NO:62) 


interferon gamma receptor 
(Ifngr) 

(SEQ ID NO:63) 


lost on transformation 1 
(L0T1) 

(SEQ ID NO:64) 


Synapsin 2 
(SEQ ID NO:65) 



Algorithm Scores 


ECS 
Algorithm 
Score 


25.7 


25.4 


24.9 


24.3 


24.3 


24.3 


CN 


CO 
CN 


CO 
CN 


22.1 


21,8 


21.5 


21.4 


Groups 1, II 
(Literature) 
and III 


25,7 


25.4 


24.9 


CO 


24.3 


14.3 


CN 




CO 
CM 


22.1 


21.8 


21.5 


11.4 


PGI II 


o 


o 


o 


o 

CN 


o 


o 


O 


CD 


O 


o 


o 


o 


o 


PGI Data, Group II: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


o 


o 


o 


O 


o 


o 


O 


CO 


O 


o 


o 


o 


o 


Human 
cells 


o 


o 


o 


CM 


o 


o 


O 


O 


o 


o 


o 


o 


o 


Human 
tissue 


o 


o 


o 


00 


o 


o 


O 


O 


o 


o 


o 


o 


o 


Unigene number 


Mouse 
Unigene 


VN 


NA 


NA 


Mm.4825 


NA 


Mm.23452 


NA 


NA 


Mm.121970 


Mm.2374 


NA 


VN 


Mm.227176 


Human 
Unigene 


NA 


NA 


NA 


Hs.2256 


NA 


Hs.137168 


NA 


NA 


Hs.104519 


Hs.9730 


NA 


NA 


Hs.239106 


Rat 
unigene 


Rn.4189 


VN 


Rn .37500 


Rn.10282 


NA 


Rn.9981 


Rn.4287 


Rn.3866 


Rn.9798 


Rn.9708 


Rn.43903 


Rn.98495 


Rn.11196 


Accession 
No. 


AA891542 


X17053 


AF030088 


L24374 


S81478 


U92564 


AA891054 


AA800197 


c^ 
<j> 

CM 
CO 
CO 


M16410 


AI639109 


AA800678 


M80804 


Gene Name 
(SEQ ID NO.) 


DNAJ (Hsp40) homolog, 
subfamily B, member 5 
(SEQ ID NO:66) 


JE product (small inducibe 
cytokine A2) 
(SEQ ID NO:67) 


Activity and 

Neurotransmitter induced 
early gene 3; (ania-3) 
(SEQ ID NO:68) ; 


matrilysin ( Matrix 
metalloproteinase 7 ) 
(SEQ ID NC):69) 


3CH134/CL100 PTPase 
(SEQ ID NO:70) 


Roaz 

(SEQIDN0:71) 


EST1 94857 
(SEQ ID NO:72) 


EST 189694 Normalized 

rat heart 

(SEQ ID NO:73) 


phospholipase D gene 2 
(SEQ ID NO:74) 


tachykinin 2 (Neurokinin B 
precursor (NKB)) 
(SEQ ID NO:75) 


rx04104s (Neurdi) 
(SEQ ID NO:76) 


EST190175 (atrophin-1 
related protein) 
(SEQ ID NO:77) 


solute carrier family 3, 
member 1 
(SEQ ID NO:78) 
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Algorithm Scores 


ECS 
Algorithm 
Score 


CM 


20.3 


20.2 


19.9 






18.8 


18.8 


18.4 


17,8 


1^ 




Groups 1, II 
(Literature) 
and III 


CO 


10.3 


CM 


19.9 


a> 

T— 




10.8 


12.8 


18.4 


17.8 




in 


PGI II 


00 


o 


CD 


o 


o 




00 


CD 


o 


o 


o 


CM 


PGI Data, Group II: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


O 


o 


o 


o 


o 


o 


o 


CD 


o 


o 


o 


CM 


Human 
cells 


O 


o 


O 


o 


o 


CD 


o 


O 


o 


o 


o 


O 


Human 
tissue 


00 


o 


CO 


o 


o 


00 


00 


O 


o 


o 


o 


O 


Unigene number 


Mouse 
Unigene 


Mm.3726 


Mm.976 


Mm.203907 


NA 


NA 


CO 
CD 

E 


Mm.101224 


NA 


Mm,6635 


NA 


Mm.7549 


NA 


Human 
Unigene 


Hs.7594 


Hs.1424 


Hs.172180 


NA 


NA 


Hs.285013 


Hs.80561 


NA 


Hs.374536 


NA 


Hs.400353 


< 


Rat 
unigene 


Rn.95055 


Rn.867 


Rn. 10835 


NA 


NA 


Rn.10123 


Rn.62694 


Rn. 40672 


Rn.147 


NA 


Rn.11183 


Rn.25 


Accession 
No. 


D 13962 


M84719 


AI237576 


S74907 


D82074 


D32209 


All 76710 


AI639256 


All 02838 


X17012 


U30186 


AA859627 


Gene Name 
(SEQ ID NO.) 


Neuron glucose 
transporter 
(SEQ ID NO:79) 


flavin-containing 
monooxygenase 1 ' 
(SEQ ID NO:80) 


SPA-1 like protein 
(SEQIDN0:81) 


PP1MM110 
(SEQ ID NO:82) 


BHF-1 

(SEQ ID NO:83) 


leucine-rich acidic nuclear 
protein 

(SEQ ID NO:84) 


Nuclear receptor 
subfamily 4. group A, 
member 3 (N0R1) 
(SEQ ID NO:85) 


rx01019s Rattus 
norvegicus cDNA clone 
(SEQ ID NO:86) 


Isovaleryl Coenzyme A 
dehydrogenase 
(SEQ ID NO:87) 


insulin-like growth factor 11 
(somatomedin A) 
(SEQ ID NO:88) 


GADD153 
(SEQ ID NO:89) 


UI-R-E0-bs-h-03-0-Ul.s1 
UI-R-EO Rattus 
non/egicus cDNA clone 
(SEQ ID NO:90) 
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Algorithm Scores 


ECS 
Algorithm 
Score 


16.8 


CD 


15.8 


15.8 


15.7 


15.7 


15,4 


15.1 


in 


14.3 




CO 


Groups 1, II 
(Literature) 
and III 


16.8 


O 


15.8 


00 

ih 


15.7 


15.7 


15.4 


ai 


in 


14.3 




CO 


PGI n 


o 


CO 
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CO 
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PGI Data, Group II: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 
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Human 
cells 
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o 
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CO 


o 
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o 


o 


Human 
tissue 


o 
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o 


o 
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o 


o 


o 


o 


Unigene number 


Mouse 
Unigene 


NA 


NA 


NA 


Mm.45396 


Mm.24662 


Mm.7320 


NA 


Mm.26378 


NA 


Mm.56904 


NA 


Mm.24193 


Human 
Unigene 


NA 


NA 


NA 


Hs,1 35578 


Hs.433503 


Hs.288261 


NA 


Hs.18791 


NA 


Hs.1 84889 


NA 


Hs.2699 


Rat 
unigene 


NA 


Rn.3018 


< 


Rn.11199 


Rn.11091 


Rn. 10636 

1 


Rn.6975 


Rn. 14758 


NA 


Rn. 10540 


Rn.6875 


Rn.7044 


Accession 
No. 


X06916 


AA866485 


S65091 


M95768 


AI228548 


U66479 


M65149 


AA892511 


All 76460 


AI230211 


AA875084 


L02896 


Gene Name 
(SEQ ID NO.) 


Calvasculin; S100 
calcium-binding protein A4 
(SEQIDN0:91) 


UI-R-A0-bd-e-03-0-UI,s1 
UI-R-AO Rattus 
non/egicus cDNA clone 
(SEQ ID NO:92) 


cyclic AMP-regulated 
phosphoprotein 
(SEQ ID NO:93) 


di-N-acetylchitobiase 
(SEQ ID NO:94) 


EST225243 (SI 00 
calcium binding protein 
AI) 

(SEQ ID NO:95) 


mothers against dpp 3 

homolog 

(SEQ ID NO:96) 


CCAAT/enhancerbinding, 
protein (C/EBP) delta 
(SEQ ID NO:97) 


EST196314 (tescalcin) 
(SEQ ID NO:98) 


EST220045 
(SEQ ID NO:99) 


Potassium voltage gated 
channel, Shal-related 
family member 2 
(SEQ ID NO: 100) 


transducln-like enhancer 
of split 1(Tle1) 
(SEQIDNO:101) 


Glypican 1 
(SEQ ID NO:102) 
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Algorithm Scores 


ECS 
Algorithm 
Score 


11.2 


10.7 


10.6 


10.4 


10.3 


o 


o 


CNi 
CJ) 


CJ) 


00 
00 


00 


lO 
00 


Groups 1, II 
(Literature) 
and III 


11.2 


10.7 


CO 


10.4 


10.3 
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CVJ 
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00 
00 


00 
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00 


PGI II 
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CD 
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PGI Data. Group II: Scores 
for gene overlaps with: 


Mouse 
+ Rat 
tissue 


o 


o 


CD 


o 


o 


o 


o 


o 


o 


o 


O 


o 


Human 
cells 


o 


o 


O 


o 
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o 


o 


o 


o 


o 


O 


o 


Human 
tissue 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


O 


o 


Unigene number 


Mouse 
Unigene 


VN 


NA 


Mm.41508 


NA 


VN 


NA 


NA 


NA 


NA 


NA 


VN 


VN 


Human 
Unigene 


NA 


NA 


Hs.296420 


NA 


NA 


NA 


NA 


NA 


VN 


VN 


NA 


NA 


Rat 
unigene 


NA 


Rn.10581 


Rn.43195 


Rn.1072 


Rn.22703 


Rn.7498 


Rn.20963 


Rn.4040 


Rn.22698 


Rn.22469 


Rn.29900 


NA 


Accession 
No. 


X51615 


All 69756 


AI639169 


AJ007016 


AA891735 


AA893980 


AI639058 


AA874943 


AA891580 


M61725 


AI009141 


S77494 


Gene Name 
(SEQ ID NO.) 


connexin protein Cx26 
(gap junction beta-2 
protein) 

(SEQIDN0:117) 


EST21 5655 (similar to 
G33) 

(SEQIDN0:118) 


n<04422s Rattus 
norvegicus cDNA clone 
(SEQ IDN0:119) 


protein tyrosine 
phosphatase (PtP4a2) 
(SEQIDNO:120) 


EST1 95538 
(SEQ IDN0:121) 


EST 197783 Normalized 
rat placenta 
(SEQ ID NO:122) 


rx04826s Rattus 
norvegicus cDNA clone 
(similar to Nedd4) 
(SEQ IDNO:123) 


UI-R-E0-ci-e-02-0-UI,s1 
UI-R-EO Rattus 
norvegicus cDNA clone 
(SEQ IDNO:124) 


EST 195383 Nomiallzed 

rat kidney 

(SEQ IDNO:125) 


Transcription factor 
UBF(SEQIDNO:126) 


EST203592 (chromobox 
homolog 1 ) 
(SEQ IDNO:127) 


Lysyl oxidase 
(SEQIDNO:128) 



CO 

00 



U 
O 

q 

CM 
00 

S 

CM 
CA 

00 



CO 



Algorithm Scores 


ECS 
Algorithm 
Score 


00 




C75 
CO 


CD 


CD 
CD 


in 

CD 


CO 
CD 


CN 
CO 


CN 

in 


CM 

in 


m 




^ ^— 

Q. CD -o 

3 0) £ 
> — 


00 




<7) 
CD 


CD 


CD 
CD 


in 

CD 


CO 
CD 


CN 
CD 


CNJ 

iri 


CN 

in 


tn 




O 
Q. 


o 


o 


O 


O 


O 


o 


O 


O 


o 


o 


o 


O 


PGI Data, Group II: Scores 
for gene overlaps with: 


CD a, 

CO CD 3 
O m 


o 


o 


O 


O 


O 


o 


O 


O 


o 


o 


o 


o 


CO w 

El 


o 


o 


O 


O 


O 


o 


O 


O 


o 


o 


o 


o 


E w 

I 


o 


o 


O 


O 


O 


o 


O 


O 


o 


o 


o 


o 


Unigene number 


(D ? 

If 


< 


< 


< 
2 


< 
2 


< 
2 


< 

2 


< 
2 


< 
2 


< 

2 


< 

2 


< 

2 




|i 


< 


< 


< 
2 


< 
2 


< 
2 


< 
2 


< 
2 


< 
2 


< 

2 


< 

2 


< 
2 


< 

2 


tr: .2^ 

c 


00 
CO 

d 


< 
2 


o 
eg 


eg 

lO 
c: 

q: 


^ 
CD 
CO 
CN 
CN 

C 


< 
2 


00 

c 

q: 


CO 
CM 

<^ 
d 


CD 
CD 

d 
IT 


CD 
CN 
CN 

d 


< 
Z 


in 
c 

a: 


Accession 
No. 


05 
CO 
CO 
CO 

o 
X 


o 

CO 

1^ 


5) 

00 

1 


CN 
CO 
CO 
CN 
C35 

i 


o 
<j> 
<j> 

CO 

X 


00 

r*- 

00 

in 
o 

LL 
< 


CN 
CN 
00 
CO 
C7) 

00 


in 

CO 

o> 

CN 
N 


in 
in 

05 
CO 
CD 

< 


00 
CD 

^ 

05 

in 
oo 

% 


in 
CO 


o 

C35 
CD 

c:5 
in 

1 


Gene Name 
(SEQ ID NO.) 


CO 

5 _ 

.Q CN 
1— V. 

c d 
V in /-y 

5^ 


CD 
CO 

cc ^ 
.9 o 

^ CO 

ii 

CJ) — 

=5« 


c 

p ^ 

^ Q-P 
03 CO 

J. c 

CD r— ^* 

liS 


CN 
CO 

in V-/ 

CD 2 

CD Q 
C3) — 

uj S2. 


CO 
CO 

d 

<^ 2 

CD Q 

o — 


CO 

<D ^ 

S5 
C - • 

§9 

E LU 

£ J2 


CN 

in 

9 5? 

i. 5 

a> o — 

ego 


^ CO 
0) ^ 

<2 rS 
J^2 

c — 

CD O 

o 


CD 

c 
o 
o 

3 2 CO 
Q 

00 O) — 


CsJ 

o c m 00 

o c C CO 

ll'ld 
o 8 

CO -o (0 ^ 


c ^ 

CO I- C35 

S.E d 
e2 oQ 

^ ^ C LU 


W CD 

c 

3 o 

< O 

CO 2 

6 o Q 

LLI lU C3) — 

55 gja 



o 
o 
CO 



o 



CO 
0 

- CO 
= Q. 
Q. i5 

O o 

CO c 

CD 
Q. 



0) 

E 
c 

0 
c 

CD 

"c 
3 



CD 



E 

CO £ 

^ cpco 
< 



— <D 

T 

Q. CO X3 

5 03 £ 

p ^ CO 



a. 



0) m 
CO (0 3 

g cr: w 



c 

CO <0 

El 



E w 

3 CO 



if 



0 
3 



C 

g 

CO 
CO 

8 



Z Q 

LU 

Oco 



o 

(DO 



O = Q 
05 7- — 

^ oO 



CO 9 



UJ 



<o 6 

O Q 
o> — 

^co 



LU CN 



CO 



X 
0) 
Q. 

E ^ 

O Tj- 

O o 'f- 

& Q.U 

txj .g — 

<D <D a 



00 
CO 



00 
CO 



CD 



2 2! 
^ CO 



CO 



CO \^ 



in 

CO 



in 

CO 



a. ^ 

03 W 



E 

CO 00 

CD "T 



CD — 



CO ^ or 



CO 



CO 



CD 

c 
g 
o 

<o? 
Q r: 

O CO 



in 

0 
C CD, 
jii CO 

<D L 
CJ) . 



cb g^O 
CD -o 52 



C3> 



csl 



c 

<D 

E 

o 



GQ 



CO 



"§ "D 

"S ^ 
CD 

o) 

if 

o ^ 

Q. O 



<D ^ 
CNJ 

E 

CD n 
EC 

19 
^ o 

^ LU 

3ei 



oo 



m 

u 
o 
q 

oo 

o 
o 

Q 



r- 
oo 
o 

e 

o 
o 

»o 

fO 
CM 



9. REFERENCES CITED 

Numerous references, including patents, patent applications and various publications, 
are cited and discussed in the description of this invention. The citation and/or discussion of 
such references is provided merely to clarify the description of the present invention and is 
5 not an admission that any such reference is "prior art" to the invention described here. All 
references cited and/or discussed in this specification (including references, e.g., to biological 
sequences or structures in the GenBank, PDB or other public databases) are incorporated 
herein by reference in their entirety and to the same extent as if each reference was 
individually incorporated by reference. 
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